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1.0 INTRODUCTION 

1.1 Identification and Significance of tbe Problem 

With the ever increasing use of space satellites, and the concomitant increase in the amount of 
man-made debris, the issue of accidental collisions with particles becomes more severe. While 
the natural micrometeoroid population (Cour-Palais, 1969) is unavoidable and assumed constant, 
continued launches of spacecraft steadily increase the amount of debris (Kessler, 1 #87 and 1990). 
The latter comprises items ranging from spent rocket cases, through defunct satellites, to the 
remains of bodies which have fragmented due to explosions or collisions. The result is a 
distribution of particles ranging in size from microns to meters. As with die m xometeoroids, 
die population of debris is greatest for the smaller particles, tending to decrease with an inverse 
power-law versus size. 

To understand and model these environments, impact damage features (/.< , craters and 
perforations) from returned spacecraft materials (e.g.. Long Duration Exposure Facility [LDEF], 
Solar Maximum Mission, etc.) must be analyzed in order to determine the impact parameters 
(e.g., particle size, particle and target material, particle shape, relative impact speed, etc.) 
associated with each feature. Such analysis requires die use of generic analytic scaling laws 
which can adequately describe the impact effects. Currently, most existing analyse scaling laws 
are litde more than curve-fits to limited data and are not based on physics, so are iot generically 
applicable over a wide range of impact parameters. Therefore, generic physics based analytic 
scaling laws still need to be developed. 

The LDEF Meteoroid and Debris Special Investigation Group (M&D SIG) is chai ered to collect 
and interpret the LDEF impact data and apply the interpretations to verify and update the 
environment and impact effects models. Since the largest sources of LDEF data are impacts into 
aluminum and FEP Teflon, the M&D SIG's efforts will be hampered without broadly applicable 
scaling laws. Consequendy, the M&D SIG is supporting these POD Associates, Inc. (POD) 
efforts to develop new generic physics-based scaling laws. The results of these efforts are being 
presented in this report 

This report summarizes die development of two physics-based scaling laws for e scribing crater 
depths and diameters caused by normal incidence impacts into aluminum and TFE Teflon. The 
report then describes equations for perforations in aluminum and TFE Teflon for normal impacts. 
Lasdy, this report also studies the effects of non-normal incidence on cratering -id perforation. 

Although FEP Teflon was flown on LDEF, TFE Teflon is studied instead sin ; the available 
experimental data (Hdrz, 1992) is for TFE Teflon and die material parameters provided in the 
literature (Harper, 1992; Moses, 1978; Dean, 1992; Rice, 1980) reveal som ; inconsistences (some 
references even fail to distinguish between TFE and F EP Teflon, which have snghdy different 
properties). However, it is expected that TFE and FEP Teflon will show similar i jsponses under 

hypervelocity impact conditions. 
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1.2 Background 


For bodies in low Earth orbit (LEO), whether orbital debris or antisatellite weapons, there is only 
a small change in speed versus altitude even out to 2000 km, and the average speed is about 7.7 
lcm/s (e g-, at 500 km). However, because different bodies are in different orbits, collisions are 
possible with impact speeds between zero and about 15.4 km/s (average about 10 km/s), 
dependent upon the angle between die velocity vectors. The two speed limits for same-altitude 
bodies are for "tail-chasing" and "counter-orbital" conditions, respectively. These possible impact 
speeds are independent of particle size and mass. 

Micrometeoroids orbit around die Sun rather than die Earth, but their orbits can intercept the 
Earth's. Since die Earth itself moves at about 30 km/s, while the micrometeoroids can have 
speeds of up to 42 km/s (for highly elliptic orbits) at die Earth’s orbital radius from die Sun, the 
results can involve collisions with maximum impact speeds from about 12 to 72 km/s (for co- 
orbiting to counter-orbiting cases, respectively). For a space platform in orbit around die Earth 
(at 7.7 km/s), the additional effect of velocity vector summation results in micrometeoroid 
collisions with speeds from almost zero to 79 km/s. In reality, very few collisions occur at die 
highest speeds, and die average impact speed is about 19 km/s, with a meaningful population 
range from a few km/s to about 25 km/s (Zook, 1990). As with die debris, the impact speeds 
are independent of micrometeoroid sizes and masses. 

While the debris is mosdy concentrated in die higher inclination orbits (above 60 degrees), die 
micrometeoroids effectively arrive from all directions. In both cases, die effect of the spacecraft 
orbital motion always results in the highest flux (impacts per presented area per time) and highest 
impact speeds on die RAM surface, for which the surface normal points along die spacecraft 
velocity vector. Consequendy, the RAM surface suffers die greatest degree of impact damage. 
The pseudo-circular orbits of die debris result in impacts on die spacecraft surfaces which are in 
die plane locally parallel to the Earth's surface below (i.e., RAM, SIDES and TRAIL), and die 
impact rate rapidly decreases as die TRAIL surface is approached, and there are very few 
collisions with die EARTH or SPACE facing surfaces. The micrometeoroids impact all surfaces, 
but for LEO, the Earth provides significant geometric shielding (a "shadow"), reducing the 
number of impacts on all faces except SPACE, but especially on the EARTH facing surface. 

Note that, due to orbital mechanics and satellite geometries, very few of die impacts involve 
"normal" collisions (impact velocity along the normal to die target surface), and thus die hits 
occur with a large spread in impact angle. It is well known that die angle of approach modifies 
die resulting penetration (Christiansen, 1992c), hence die effect needs to be properly understood 
in order to correctly interpret the environment and die consequential impact effects. Likewise, 
die effects of density of both die impacting particle and the target must be known. Presentiy, it 
is assumed that most of the micrometeoroids have a density of about 0.5 g/cm 3 (Cour-Palais, 
1969), while the debris value is about 4.7 g/cm 3 for small particles (less than 1 cm), and 
decreases with increase in particle size (Kessler, 1987) (since most large pieces are not chunky 
pseudo-spheres, but rather odd-shaped items which, on average behave as if partiy porous). 
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Furthermore, the effects of impactor and target material properties ( i.e ., yield, f racture strength, 
melt energy, etc.) must be understood. 

1J Damage Modes 

The impact damage ranges from simple pitting, erosion and cratering for impact into plastically 
yielding materials, through conchoidal and star cracking for brittle targe ;, to complete 
perforation, large-scale spallation/fragmentation, and material melting and vaporisation. The data 
from LDEF indicates additional effects, such as delaminations of multilayered materials, and the 
generation of rings of ejected material, and/or permanently deformed material. These effects are 
in addition to thermal cycling, UV (ultraviolet light) embrittlement of plastics, and AO (atomic 
oxygen) erosion. All these effects can be synergistic. Impacts can modify material properties 
or morphologies, thereby altering the subsequent responses to UV, AO or thermal cycling. 
Likewise, thermal cycling, UV, and AO can alter properties such that the issultant impact 
cratering, cracking and perforations are modified. These various effects, and the ' synergies, are 
all being studied. 

1.4 Issues 

To correctly interpret the space environment, the resulting modes of impact c mage, and the 
potential methods that may be employable to mitigate the effects, it is necessity to property 
understand the ’Yules" of impact damage. Unfortunately, existing experimental facilities cannot 
replicate the complete range of conditions. Gas guns are limited to maximum speeds of about 
8 km/s, electric guns have limitations due to sabot requirements. Van de Graaff accelerators can 
only be used with small charged particles, etc.. Further, many techniques have been developed 
primarily to accelerate flat plates to study planar, one dimensional (1-D) strain, impact effects. 
However, these approaches are not easily adapted to the problem of throwing :hunky bodies. 
Consequently, only partial experimental testing can be done at present. In particular, there are 
no satisfactory means of replicating the high speed (above IS km/s) impacts of the larger particles 
(above 0.1 mm diameter), of either metallic or (especially) non-metal lie nat jre, in a well- 
controlled manner that also allows for non-normal impacts. 

As a result, reliance must be put on either extrapolation of (limited) experimental data obtained 
at lower speeds, or with non-typical materials, or put on computational approaches. Although 
computational capabilities have increased dramatically during die last decade {more rigorous 
physical modeling, and much faster and more capable computers), die fact remains that properly 
calculating impact problems is time-consuming and cosdy. Furthermore, the accuracy of the 
calculation is limited by adequacy of known material properties and behavic . such as high 
pressure EOS (equation of state), rate-dependent yielding and fracturing, etc Much of the 
experimental data needed is frequentiy obtained from 1-D experiments (plate s ap). However, 
while this approach eases die interpretation of material responses, it does not necessarily provide 
all the required data. For example, 1-D strain is an irrotational flow problem, whereas cratering 
and penetration of one finite body into another clearly involves material rotational flow. Thus 
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computational techniques are usually best reserved for analysis of specific cases, and for 
identifying sensitivities to parameter uncertainties. 

As a result, a need still exists for analytic scaling laws which can adequately describe the impact 
effects. Done properly, such laws are very useful for determining the "ball-park" of probable 
responses, which in turn allows attention to be applied to specific cases via the computer codes. 
These laws need not be perfect (10 percent accuracy is sufficient for many purposes), but they 
do need to be rational and based on physics, such that extrapolations to experimentally 
inaccessible conditions are credible and well-founded. 

Unfortunately, many existing scaling laws are little more than curve-fits to limited data, and are 
usually derived by investigators who invariably fit only their own data. Accordingly, there are 
several such laws, and simple inspection reveals glaring discrepancies. For example, these laws 
rarely agree on such issues as die proper power-index to be applied to impact velocities or 
densities. Further, while some authors invoke material "hardness” numbers, others use yield 
strengths. However, invariably such parameters are applied to only the target materials, but not 
tiie impactors, despite the problem being one that deals with both bodies. Some laws fail to 
acknowledge the well-known "supralinearity" observed in cratering, wherein the crater size 
increases faster than the impactor size (all other factors remaining constant). Unfortunately, those 
laws that do incorporate this effect are themselves in violation of "dimensional analysis". 
Frequently, these formulations involve numerical "constants” that cannot be (!), since to adhere 
to physical principles they must have an inverse dimension to that "left hanging" within the 
equation. Further, while many such laws define an upper limit to validity (e g., in velocity), 
rarely is the corresponding lower limit given, despite the fact that the formulations are self- 
evidently incorrect for low velocity (or density, etc.). 

Many other problems abound. For example, if material yield strength is invoked, the value used 
is tire static one, whereas it is well known that many materials display considerable strain and 
strain-rate sensitivities, and can also suffer from thermal softening effects. Most investigators 
assume that the "cosine law” applies to impacts (i.e., the penetration is a function of the normal 
component of the impact velocity). However, while this effect is observed for the lower impact 
speeds (indeed, ricochets can occur for angles of incidence greater than about 60 degrees from 
the surface normal), it is not necessarily true at the higher speeds. This arises because, at the 
lower velocities, the impact phenomenon is dominated by momentum, but as the velocity 
increases so the process becomes dominated by energy. Thus an impact by a small particle at 
very high velocity will behave very much like a surface explosion. Such explosions do not obey 
the cosine law. Also, few investigators explain what happened to the remaining component of 
the incident energy and momentum. 

Impact events fall into three basic classes: flat plate (1-D) planar impacts, chunky (pseudo- 
spherical) impacts, and long rod impacts. For the smaller man-made space debris and 
micrometeoroids, it is reasonable to assume the particles are chunky. However, for the larger 
(and far more deadly) particles it is unlikely that they are spherical. Rather, the pieces will be 
parts of thin plates, or antenna booms (etc.). These impactors will produce results more typical 
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of those for flat plates or long-rod penetrators. Surprisingly, although many in estigators have 
studied these responses independently, there have been few attempts to many the results. Thus 
no analytic law presently exists that allows description of die gradual trend from flat plate 
throu gh spherical shape to long rod impact. Consideration of this issue will be useful in its own 
right, since it allows a simple appreciation of die consequences of impactor ch inges in shape. 
Further, this approach immediately reveals the source of confusion in die multiple power-laws 
that exist. 

1.5 Technical Objectives 

The technical objective for this work is to develop a set of generic scaling laws hat adequately 
describe die physics of impact events. This will include cratering, cracking and spalling, 
perforations, and melting/vaporization phenomena. The laws will take into account material 
properties, sizes, shapes, velocities and angles of incidence, and will logically extrapolate to both 
high velocity conditions presendy inaccessible, and to (known) low velocity conditions. 

Specifically, the objectives are to: 

(1) Study the responses of aluminum and Teflon (specifically T) E) for normal 
impacts by chunky bodies, developing a generic scaling law that describes 
cratering through perforation, and 

(2) Expand these studies to include non-normal impacts. 

This report specifically studies the effects of both normal and oblique impacts int< aluminum and 
Teflon targets, and also gives some data for other materials. 

1.6 Technical Approach 

POD has tackled die issue of impacts into aluminum and TFE Teflon vis three approaches: 

(1) Application of physical logic in an attempt to determine which pa, imeters should 
apply to the problem, 

(2) Application of hydrodynamic code calculations, using the CTH cc Je (Bell, 1991) 
from SNLA (Sandia National Laboratory, Albuquerque), to map o>it the predicted 
responses as functions of input parameters, and 

(3) Study of existing scaling laws and comparisons with both experimental data and 
the results of CTH calculations to determine which (if any) such l; w fits the data. 

This report first addresses the results of each of these approaches separately, the combines the 
results into sets of conclusions for this study and recommendations for future ei iorts. 
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2.0 PHYSICAL LOGIC 


2.1 Hypeivelocity Impact 

Exactly what constitutes a hypervelocity condition is not well defined. For some investigators 
the condition applies for all impact speeds above 5 km/s, while for others it occurs when u, > 
Cot where Uo is the impact speed and Co is the material (usually the target) stress wave speed 
(whether this speed should be the Hugoniot shock, the low-pressure bulk, the extensional, the 
longitudinal [dilatational], or die shear velocity is rarely made clear: most researchers reference 
die low pressure bulk sound speed). Another definition would be whenever gross plastic flow 
occurs in a rotational flow pattern, and the stress approximates the Bernoulli law for fluids. Note 
that this definition is obviously material-specific, being easily reached for soft ductile metals but 
applying only for higher impact speeds for the less ductile, high yield strength, ceramics and 
glasses. Also note that this definition never applies to die well known 1-D plate slap condition, 
since die true 1-D case does not involve rotational flow. POD is of the opinion that the latter 
definition is the more logical one, since the impact community refers to hypervelocity conditions 
associated with impactors of finite lateral dimensions, which invalidate 1-D stress logic (and 
make die analysis much more difficult). 

The reference to the Bernoulli stress state and associated material flow needs to be qualified. As 
normally used in die impact community (and in this report) it is assumed that the materials are 
basically incompressible. Strictiy speaking, this approach is merely an approximation which 
applies once the stresses decrease to much lower values than die initial impact stresses. Section 
3.0, which discusses the CTH calculations, elaborates on this subject. 

22 Basic Phenomenology 

We first describe the fundamental physics of impact cratering in order to establish the overall 
phenomenology. For convenience we chose our "chunky projectile" to be a cylinder striking the 
target in the manner of a 1-D plate slap, with the end face of the cylinder parallel to the target 
surface. The cylinder has length L and diameter dp, with L = d,,, and the impact speed is ^ and 
is normal to die target surface. 

Immediately upon impact two shock waves are generated, one propagating into the target and one 
propagating back into the projectile. The initial stress, o, is given by the 1-D condition 

o = Z p Z l u«/(Z p + Z t ) 0) 

where Z is the acoustic impedance given by 

Z = p(c 0 + su)u 

where p is material density, c 0 is the low-stress bulk sound speed, u is the induced particle 
motion change in speed, and the term s is material-specific and related to the material Gnmeisen 
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parameter (see Section 2.3.5). Note that for low impact speeds the stress is proportional to this 
speed, while for high impact speed the stress becomes proportional to the square of the speed. 
This 1-D stress can be very large. For a symmetric Al/Al impact, with p = 2.7 g/cm 3 , c 0 = 5.1 
km/s, s = 1.4 and an impact velocity (Uq) of 10 km/s, the stress is 1.63 Mbars ( 63 GPa). This 
stress is grossly in excess of the plastic yield strength (Y) of aluminum, which is about 2.7 kbars 
(0.27 GPa) for A1 6061-T6 (Shackelford and Alexander, 1992). 

Release waves are generated at the edge of die projectile (free surface relief) an at die edge of 
the impact region within the target. The projectile release waves cause a drop in the stresses 
within the projectile, such that after a time t = r/c,, (with r p die projectile radius, and c p die radial 
release wave speed) die 1-D shock state reduces to a value closer to die BemoJli stress. The 
stress region in the target direcdy under the impact footprint itself propagates stress laterally, 
while die corresponding release waves drop die stress under die impactor towards the Bernoulli 
state, in the time t = r/c* where c, is the release wave speed in die target. Initially, if die 1-D 
stress is o 0 the target lateral stress is (a 0 - Y). The downstream particle motion under die 
impactor causes a shear wave to propagate laterally into the target The corresponding maximum 
propagated shear stress is limited to the yield value. Additionally, the compressive lateral stress 
itself propagates a stress of value 0.5(q,-Y), since die boundary under the edge of the impactor 
is one of like-to-like material. This latter stress itself decreases with distance owing to the radial 
divergence. This stress propagates as a longitudinal wave laterally and is ini dally far more 
imp ortant than the shear wave. The radially moving longitudinal wave itself has a lateral stress 
(that is thus perpendicular to the target free surface outboard of die impactor) or (q/2 - 3Y/2). 
Upon reaching die target free surface this stress causes outer motion (upstream, and generates 
a release wave traveling back into the target Clearly, when o 0 » Y the material behaves 
essentially as a fluid with stresses approximately equal in all directions (i. s., we ave a pressure 
rather than stresses). 

Thus after the release waves cross the impact axis the 1-D stress rapidly vanishes and the target 
material is set into a rotational flow, in a similar manner to the flow of a fluid around a moving 
body. Material on the target surface outboard of the impactor flows in a combined radial- 
outward/backstream-upward manner. This is the flow associated with the formation of jets of 
material and the onset of the formation of the lips observed around erasers in ductile targets. 
Once this rotational flow is established the stress state on die impact axis becomes the Bernoulli 
one, given by 


o = 0.5 pu J (3) 

where u is the local particle speed. For the symmetric Al/Al impact this stress s given by 

o = P p UoV8 * ft Uo J /8 (4) 

which is 337.5 kbars (33.75 GPa) for an impact speed of 10 km/s. Thus the Bernoulli stress is 
much lower than the initial 1-D shock stress, but still well in excess of the aluminum yield 
strength. This Bernoulli stress is a transient state, only, for the case of a chunk, projectile and 
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itself rapidly decays unless continued momentum arrives at die projectile/target interface as 
occurs for a long rod penetrator (LRP). 

23 Cratering Behavior 

It is well known that, when a projectile impacts an infinitely thick target at high velocities, a 
crater is formed in the target. Cratering behavior for hypervelocity impacts has been studied 
emp irically for over 40 years. However, many issues have remained unresolved, including 
whether foe craters are fundamentally a function of impactor kinetic energy or momentum. 
Likewise, foe appropriate material parameters governing foe responses have beat uncertain. 

23.1 Melt Craters 

Some investigators assume that craters in targets are associated with foe development of impact- 
induced melting. However, this implies that foe cratering phenomenon ceases once foe stresses 
drop below those necessary to promote melting caused by foe excess entropy trapping of shock 
waves. This would cause predictions of craters to be noticeably smaller than those observed, 
especially for very ductile targets. The true formation of craters must therefore be related to 
plastic yield flow. This makes sense since plastic flow is merely a fluid flow with off-set stresses 
(i.e., stress deviators exist in foe solid material). Melt and vaporization flow is merely plastic 
flow without stress deviators (we here ignore deviators due to material viscosity, and also ignore 
surface tension phenomena. Note: almost all hydrodynamics codes also ignore these two 
phenomena. Hydrocodes include artificial viscosity, but this is done only to control numerical 
instabilities that frequently occur when computing strong shocks). If foe target cratering involved 
pure melt (i.e., pure fluid) conditions another problem could also arise, namely foe self-healing 
of foe crater due to foe internal pressure within foe fluid and foe effect of surface tension. This 
latter behavior is typical of impacts into true fluids (e.g., into water), where it is well known that 
a transient crater is formed that disappears later due to foe inflow of foe material, leaving no 
permanent crater. 

Clearly, a "melt" crater is foe minimum possible crater, and can be expected to dominate foe 
cratering process only for cases where foe yield stress of foe target is very large and foe crater 
surface is rapidly decelerated once foe propagated stresses drop below those that induce melting 
upon release. Thus foe response is also sensitive to foe rate of decrease in foe yield strength 
versus temperature (thermal softening). Materials which soften rapidly with increase in 
temperature (e.g., ductile metals) will effectively extend foe apparent "melt" region, whereas 
materials which soften significantly only as melt is approached (e.g., ceramics) will restrict the 
apparent melt region. 

For reference, it should be noted that for foe symmetric aluminum/aluminum impact incipient 
melting occurs for an impact speed of about 5.6 km/s (1-D stress of 0.65 Mbar) while complete 
melting occurs at about 6.6 km/s (1-D stress of 0.9 Mbar). Incipient vaporization occurs for an 
impact speed of about 10.2 km/s (1-D stress of 1.67 Mbar). Thus no melt crater can be formed 
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for impact speeds of less than 5.6 km/s for this case, yet significant cratering h experimentally 
observed in aluminum for impact speeds less than this value. 

2.3.2 Initial Approach - The "Eneigy Equation" 

One well known approach to cratering is the energy equation, which equates th kinetic energy 
of tiie impactor with the energy necessary to excavate a hemispherical crater. Titus we obtain 

(jc/ 12) ft, dpV = (Jt/12) ft d, 3 E, (5) 

where p is the density of projectile or target, d is the diameter of projectile or rater, ft is the 
impact speed, and E, is a limiting energy-per-mass to create the crater wi tin the target. 
Rewriting, we obtain 

d^dp - (ft/ft) 1 * ft M / (E.T (6) 

Thus this equation gives the familiar 2/3rd law for velocity, and a l/3rd 1. v for density. 
However, there are obvious problems with this approach. Energy absorption (f ^r mass) in the 
target is presumed to be a constant within the crater, which is clearly wrong s ice the energy 
density varies with the stress levels. Further, the logic ignores energy absorpdon within the 
impacting projectile, and also ignores the elastic energy stored in the stress wave propagating 
beyond the crater region. Thus the energy used to excavate the crater is always less than the 
original impactor kinetic energy, but the fraction varies according to the relative strengths, 
moduli, thermal properties of the projectile and target, and the relative mass of the crater versus 
the projectile. The use of a limiting energy density is really the same as defining a limiting 
stress, since energy-per-mass is equivalent to stress-squared divided by twice a modulus and 
divided by the density. 

Thus E, = 1/2 ft 2 /(ft M) (7) 

where ft is the limiting stress, and M is a modulus of the target Despite thes problems this 
energy equation does give a reasonable account of cratering provided an empiriejd multiplier is 
used. For example, using the incipient melt enthalpy of aluminum of 150 cal/g (6 28 x 10 9 erg/g) 
and an impact speed of 10 km/s. Equation 6 predicts d/dp = 5.42 which is cli ;e to observed 
values. 

133 More General Approach 

Assume that upon impact a pseudo-spherical compressive wave is propagated into the target 
centered on the impact point. Release waves from the target surface outboard of the projectile 
will cause the compression to be followed by a tensile wave. Release waves from within the 
impactor do not cause tensile states by themselves since there is no inter-body strength between 
the impactor and the target surface. Hence any attempt to propagate a tensile wave would merely 
lift the impactor away from the target. 


9 


For hypervelocity impact the initial stresses and shock-induced heating will cause gross plastic 
flow in the target. If the heating is sufficient to induce a state above melting (including 
vaporization) die material will flow as a simple fluid. Once die stresses drop into the solid 
response region the flow will be standard plastic yielding. TTiere are three cases of such yielding: 

(a) during die passage of die compressive pulse, 

(b) during die passage of the trailing rarefaction pulse, and 

(c) after the main pulse has passed a given point owing to the induced radial 
divergence motion which leaves die material in a state of hoop strain and thus 
hoop tension. 

If the wave motion were truly 1-D, the radial (i.e., longitudinal) stress would reduce to zero while 
the lateral stresses would attain a permanent (residual) stress due to plastic hysteresis. For a 
radial stress greater than 2Y(l-u)/(l-2u) (i.e., twice the Hugoniot Elastic Limit, see Section 2.3 .7) 
this permanent stress is the yield value and is compressive (u is Poisson's ratio for the target). 
For a rad ial stress between the one quoted above and Y(l-u)/(l-2u) the permanent stress lies 
between zero and Y. For lower radial stresses the behavior is elastic only and no permanent 
stress occurs. 

However, die response is not 1-D and involves divergence. Thus die true lateral stresses (i.e., 
die hoop stresses) have die tensile hoop stress superimposed, giving a net tensile stress towards 
die impact point. Only for radial distances far from die impact point will these stresses become 
compressive as for the 1-D case, since the hoop tensions will become very small. These effects 
are seen in die CTH calculations discussed in Section 3.0. 

Assume that the latest state (c) defines the final crater surface. The divergent flow-induced hoop 
strain is given by Ar/r where A r is the induced radial motion at a distance r from the impact 
point The corresponding hoop stress then becomes 

ae = p, c, 2 Ur dt/r * o, dt c/r ( 8 ) 

where o e is die tensile hoop stress, a t is die radial stress for die pulse, p, is die target density, 
Cj is die target sound speed, u, is the induced particle speed and dt is die pulse time-width. Note 
that the qu an tity c^dt represents the outward momentum per unit area at position r. 

For a spherical projectile die mean time for a full release wave reverberation across the projectile 
is given by dt = 8r,/3c p where r p is the projectile radius and c p is the projectile sound speed, and 
we assume that this time defines the effective width of the radial pulse. 


Assume that die radial stress obeys die momentum rule 
a r r 2 = o 0 r p 2 


( 9 ) 
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where q, is the stress generated at the interface between the projectile and the arget. Further, 


««nnu» that this latter stress is given by the Bernoulli law such that 

O 0 = 0.5 p Uo J = 0.5 ft u, 2 = 0.5 p p (Uo - u,) 2 (10) 

where Uq is die impact speed and u, is the induced target particle speed. 

Thus p - P p ft /(P p 1/J + ft 1/J ) 2 - (11) 

Assembling die terms we obtain 

Ob = (8/3) Oq (r p 3 /r 3 )(c/c p ) (12) 

Now set this hoop stress to the yield strength of die target, Y t , and set the com ponding radial 
distance from the impact point as the crater radius, r e , to finally obtain 

(r/r p ) 3 = (4/3)( Pp /ft)( P /Y t )(c/c p )u 0 2 /(l + ( Pp /ft) ,/2 ) 2 (13) 

Thus 

d^dp « rjr p = 1.10064( Pp /p 1 ) 1/3 (p/Y t ) ,/3 (c f /c p ) 1/3 u, w /(l + (p (14) 
Several interesting points should be noted about this equation: 


(1) There is no numerical multiplier which did not come frorr the a ilysis. 

(2) The equation "looks like" the well known "energy equation", bu was based on 
momentum. 

(3) The velocity index is the familiar 2/3rd value. 

(4) The density index is not a constant value, but depends on the rat ) of densities. 

(5) The crater size depends on die ratio of target density to target yield strength. This 
term in the equation has die same dimensions as an inverse sound speed to the 
2/3rd power, but is not a material sound speed. Nevertheless, the :erm serves the 
same function as a sound speed, in a similar manner to that invoked by others, 
who frequendy use terms of the form (u</c) in their scaling laws. The yield 
strength is that of the target only, and does not need to be referenced to any other 
material. Further, note that the yield value used is die low-stress static one. The 
rationale for this is that the model describes the terminal conditions of crater 
formation. Factors such as strain-rate hardening and thermal softening are 
subsumed into the physics of die more highly stressed regions of the crater region. 
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and are treated as energy absorbing mechanisms within this region. (Note: our 
logic does not invoke energy as a parameter). 

(6) The equation is sensitive to the ratio of target and projectile sound speeds. 

(7) The equation obeys die rules of dimensional analysis, since the grouping is 
dimensionless. 

The density function can be replaced by a simple ratio if we use the equivalence 

((VP«) N = (P/Pr /(l + (Pp/PO'T 3 05) 

Doing this we immediately find that the index N varies. For a small ratio of projectile density 
to target density N tends to 1/3. As die density ratio increases so N increases. When the ratio 
is unity N is indeterminate. For a ratio larger than unity N becomes negative, being large near 
die unity ratio, decreasing towards zero for very large ratios. This "odd" behavior may be the 
"excuse" for the varying indexes found in die literature. 

Note that most experimental conditions prevent die choice of individual material parameters. To 
change one property (e g., density) invariably involves changing other properties (e.g., yield 
strength, sound speed, etc.). Thus experiments invoke "clusters" of properties. For this reason 
dependencies on such factors as material sound speeds are difficult to identify. 

For die well-known symmetric impact of aluminum into aluminum, with density of 2.7 g/cm 3 , 
sound speed of 5.1 km/s, yield strength of 2.7 kbar (i.e., A1 6061 T6) and Uo * 10 km/s, we 
obtain djd^ = 6.934. This is near the observed values. However, this is probably a coincidence! 
This case is also one that "removes" many dependencies owing to die unity value of the groups 
involving density and sound speed. 

POD believes that Equation 14 contains all die pertinent material parameters for hypervelocity 
impacts into highly ductile targets. Surface fracture craters are not presendy included in this 
model owing to the assumption that highly ductile materials consume significant energy during 
die strain-to-fracture. This effect causes a rapid decrease in peak stress and thus prevents fracture 
(i.e., it results in large fracture toughness; see Section 2.3.10). However, die equation has several 
pitfalls and die power indexes are not necessarily correct. 

23.4 Caveats 

The caveats include: 

(1) The assumption of a truly hemispherical stress wave is not correct. This point is readily 
demonstrated by considering the impact of a chunky cylinder as described before. The 
expanding compressive stress front only asymptotes towards a spherical wave at radial 
distances which are large compared to the impacting projectile radius. Thus for crater 
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radii which are not large compared to impactor radius die wave front s unlikely to be 
spherical Similar arguments apply to the target surface-generated release waves which 
cause tensile waves to follow the compressive pulse. These waves are centered outside 
of die impactor, and only coalesce into a single pseudo-spherical wav after travelling 
several impactor radii. 

(2) Eq uati on 9 describing stress versus distance presupposes a constant pul }-width. This is 
only true for a non-dispersive condition (e.g., a linear-elastic response. In reality, the 
dispersion associated with hydrodynamic wave propagation always ca ses the pulse to 
widen. This has the effect of modifying die power index for stress ve ;us distance. 

Thus we have ^ r 2 dr r = ^ r p J dr p 0^) 

where dr is the pulse-width at either the radius r from the impact point, c at the impactor. 


Dispersion gives dr r > dr p . Thus a more correct description would be 

q i* 1 = a*, r p N with N >2 0^) 

This will modify Equation 13 to give 

( r< /r p ) N+1 = (4/3) (p/pi) (p/Y,) (c/c p ) Uo 2 / (1+ (Pp/ft) ,/2 ) 2 ( l8 ) 

then we obtain the following equation 

( r ^4/3) w *'wp^v/v/ m wc/^v w *'V(i-K^prr 0 (i») 

Thus we finally obtain 

(dyd^ttC^Pp/p^ + (P./POT™ ( 20 ) 


Note that this modifies the various power indexes. The discussion abo e also defined N 
> 2. As an example, if N = 2.5 

(d^<V-1.0»57(p/pi‘“ , (p/Y 1 )““(c t / S )““«. ,J,, ‘/(l + (P/P ) m r n4 <2» 

Thus a proper description that includes the dispersive nature of the prop ^gating wave will 
give a more complete equation for cratering in ductile targets. Note that N > 2 produces 
a power index for impact velocity which is somewhat smaller than 2/3. An index of close 
to 0.58 has been previously suggested by Holsapple (1987). For this eason we choose 
to "lock" our equation as given in Equation 21. 

Substituting values for A1 6061-T6 for both the projectile and the tar et, and assuming 
an impact speed of 10 km/s, we obtain djd, = 5.26 which is also clo e to observations 
and to the results from the CTH hydrodynamic code calculations (see below). 
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(3) The logic presupposed that the active driver-stress ceased after the time t = 8r,/3c p , which 
is one complete reverberation through die projectile, both axially and laterally. This does 
not account for die entire projectile momentum. However, it probably does describe die 
conditions for die diameter of a crater since the surface of die target is die region most 
affected by the Bernoulli flow. The late-time remaining momentum does not contribute 
to the lateral radially propagating stress, rather it contributes to die axially propagating 
stress. This is because once die full diameter of die projectile has penetrated die target 
surface, diere is no further lateral "push" to the target. This is one reason why craters are 
not truly hemispherical. The same logic is expected to apply to die crater diameter 
produced by an LRP. For an LRP the initial stress conditions are identical to those of a 
chunky projectile. Only after the time 8r/3c p does die target "know" that "extra matter" 
exists in the rod. Thus beyond this time the Bernoulli stress at die head of the rod 
remains. Use of the Bernoulli law leads to die well known penetration law for LRP 
(Eichelberger, 1936). Thus rewriting Equation 10 gives 

u/(uo-u) = (Pp/pJ ,/2 (22) 

and if we assume the time of action is 

t = L/(Uo - u) (23) 

where L is the length of the rod being "consumed" with a speed (u«, - u), then the 
penetration (P) is given by 

P ■ u t « (pp/p,) ,/2 L (24) 

which is the equation describing the penetration of either a jet or an LRP into a target 
(ignoring die initial and final states of penetration, and material strengths). 

(4) The equation includes the two stress wave speeds c p and c,. However, these speeds are 
not simply the low-stress values, but are averages over die entire stress histories. A more 
accurate integral needs to be performed in order to obtain die appropriate values. 

2.3.5 Crater Diameter Venus Penetration Depth 

The above logic is thought to apply to die development of crater diameters (d c ). There is no a 
priori reason to expect the rules for penetration depth (P) to be the same. Part of the reason is 
that, whereas it could be justified to include only part of the projectile momentum for the lateral 
motion diere is no obvious logic to use anything other than die total momentum for describing 
the penetration depth, particularly if hoop-stress logic is employed. This immediately leads to 
a different answer. 
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Again assuming a pseudo-hemispherical crater, with equal local momentum p r area over the 
entire hemispherical surface regardless of angle relative to the impact axis, the axially-directed 
momentum per area at a distance r (on axis) from die impact site is given by 

I = (4n/3) p p r p 3 Uo/nr 2 ( 25 ) 

Using the hoop stress logic then leads to 

P/dp = (1/2) (4/3) tn (pp/p,) 1 * (P/Yf (quo) 10 (26) 

Thus this equation gives a direct l/3rd law for densities, but also gives a l/3r« law for impact 
speed. A modification in this law can be invoked by noting that die sound s eed term, c„ is 
really stress dependent, and should be c, = c* + su where u is the induced loc l particle speed 
in the target. 

Assuming the Bernoulli state to apply, we have u = Uo /(I + (p/p p ) ,/2 ) Thus 

P/dp=( 1 /2)(4/3) 1/3 (Pp/p J I/3 (p /Y,) 1/3 { (Cof^sUo/( 1 +(ft/f^,) l/l ))Uo} 1/3 ( 27 > 

For an Al/Al (6061-T6) impact at 10 km/s this equation predicts P/dp = 5.86 w ich is too large 
by a factor of almost two compared to both experimental data and the results front the CTH code 
calculations (see later). We have not yet found the reason for this factor. To make the equation 
more in keeping with experimental data we insert a factor of 1/2. However, th s factor is made 

peraument 

Next, we take note of material strengths (see below) by introducing a lower limit velocity term, 
u^n, (target, critical speed). Equation 27 then becomes 

P/dp=(l/4)(4/3) 1 °(Pp/pJ ,/3 (p/Y f ) l/3 {(c <t +s(u,-u tiak )/(l+(p f /pp) ,/I ))(u, u**)} (28) 

Note that Equation 28 gradually transforms from a l/3rd law for velocity into a /3rd law at high 
speeds, becoming 

P/dp =o(l/4) (4/3)‘° (Pp/ft ) ,/3 (p/y t ) ta (uo) M (s /(I + (p/p P ) ,/2 )) 1/3 ( 29 ) 

This transition depends on the ratio s/c^ with a low ratio causing the response to remain in the 
l/3rd law regime even for moderately high impact speeds, while a high ratio transforms into the 
2/3rd law at lower speeds. Theoretically, the value for s = (1 + T)/2 where T is the Gruneisen 
parameter. Materials with very low F thus have low values of s. Such materials are frequently 
porous and highly energy absorbing (r.e., dispersive). Most solids have high values of F (about 
2.0) and absorb energy less (although shock entropy-trapping still causes a degree of absorption). 
We thus expect very energy-dispersive targets to display the l/3rd law at moderate impact speeds, 
while most solids will rapidly display the 2/3rd law. These expectations are in accordance with 
the conclusions of Holsapple (1987) who states that very energy-absorbing syst ms tend to obey 
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die momentum rule for penetration depth, while truly non energy absorbing targets tend to obey 
die energy rule. The meaning of energy "absorption" is that part of die projectile kinetic energy 
is irreversibly converted into heat and/or plastic deformation, which always occurs with porous 
materials and also occurs to a lesser extent with strong shocking in solids. Note that for shocks 
sufficiently strong to cause complete vaporization die heating is recoverable, since the energy 
eventually is converted into kinetic energy of gas motion. 

Note that if the local momentum is not uniform versus angle from die axis this merely changes 
die leading numerical multiplier. For example, if momentum obeys the rule 1(9) = I(O)cos0 
then the term 4/3 becomes 2. 

As with Equation 21 which described the crater diameter, POD believes that Equation 28 contains 
all die relevant material parameters for describing crater depth, but that die various power indexes 
are possibly incorrect. 

The u^ term is not included in Equation 21 for crater diameters. The logic for this is that the 
near-free surface target response involves strong release motion. This means that the distinction 
between elastic-plastic response and the true Bernoulli response is less apparent, so die transition 
is more gradual. 

2.3.6 Projectile Strength 

The previous discussions have referenced only the target strength, but not die projectile strength. 
During die initial pseudo 1-D impact the stress states in both die projectile and the target are 
functions of their respective Hugoniots. At low stresses material strength results in an elastic 
region up to die Hugoniot Elastic Limit (HEL) and an elastic wave only is propagated. For 
stresses above HEL but below a stress d (defined below) a two-wave structure develops, 
consisting of an elastic wave precursor followed by a slower plastic wave. In this stress regime 
die precursor can pre-condition die material for later response to die plastic wave. For stresses 
above d the plastic wave overtakes die elastic wave, resulting in only one shock wave equal to 
die plastic wave. The stress d is defined as that point on the Hugoniot given by the condition 
that die Rayleigh line from the HEL up to die Hugoniot be an extension of die path from zero 
stress to die HEL. Generally, both the HEL and d are much smaller than the impact stresses for 
hypervelocity conditions. 

Upon passage of the shock wave through the projectile to its rear surface, stress release will occur 
via a rapid elastic decrease and a slower plastic decrease (unless the liquid or vapor state is 
induced, for which both yield strength and shear modulus vanish), and these release waves 
propagate back to the projectile/target interface. However, from this time on the exact state of 
the projectile has minimal effect on the subsequent cratering process, while the strength of the 
target continues to influence the cratering process. This is die primary reason why only the target 
strength is referenced. 
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23.7 Strength Effects at Low Impact Speeds 

However, these arguments are only valid provided the impact stresses are much higher than the 
HEL and o values. For low speed impact die material strengths become important, ultimately 
do minatin g the response. Clearly, if die impact produces stresses below die HE] of die projectile 
no permanent deformation of the latter occurs. Likewise, if die stress is below the HEL of the 
target no cratering occurs (not even simple plastic indentation). The situation is similar to die 
correction terms used by Tate (1967) for LRP, where die Bernoulli equation is modified by 
adding a strength term to both die projectile and the target, as in 

o = 0.5 p^, (uo - u) J + Y„ = 0.5 ft u* + Y t (30) 

Note that this correction is itself wrong since die equation does not allow for u solution at low 
values of velocity. A more correct version would treat die strength terms in ae form HEL = 
p m^ and Y = O.Spu^ 2 , with u^ setting a velocity limit for each material, his form of the 
equation does allow for solutions as impact speed drops to zero. 

Specifically, if the Bernoulli stress limit is not achieved there will be no "fluid f! w" induced and 
die response will become elastic-plastic only. If the HEL is not achieved the n ,-ponse is purely 
elastic. The critical velocities for Bernoulli yielding are given by 

for projectile: u, OT = - (2Y,/p 1 ) 1 ' 1 (l-Kiyp,)“) (31) 

for target: n. - (2Y/p f> ) , "( l+fp/ft) 1 ' 5 ) - (2Y/p^l-t<A/fc)“) ( 32 ) 

For an Al/Al impact of A1 6061-T6, die critical Bernoulli speed is 0.8 *4 km/s. The 
corresponding speed limit for the HEL is much lower, being 0.039 km/s. Thu between these 
two impact speeds the response is elastic-plastic, and no Bernoulli flow occurs 

For TFE Teflon, Y, is about 300 bars and p, is 2.17 g/cm 3 , while for uluminm 6061-T6 the 
corresponding values are 2.7 kbars and 2.7 g/cm 3 , respectively. Thus for an aluminum projectile 
and Teflon target, u^ is 0.315 km/s, while u^ is 0.946 km/s. 

As an example of these effects, consider an impact between an A1 6061-T6 projectile and an 
alumina target. For alumina the density is 4.0 g/cm 3 and die yield strength is about 60 kbars. 
We find the critical Bernoulli impact speeds are 3.84 km/s for flow in the alumina and 0.815 
km/s for flow in the Al. Thus an impact speed of less than 3.84 km/s will not induce Bernoulli 
flow-style hypervelocity craters in the alumina target, but will cause flow in ti e Al projectile, 
murin g a crater which is shallow relative to its diameter. The reverse occurs if projectile and 
target are exchanged (i.e., the projectile does not suffer flow but die target does, giving a crater 
which is deep relative to its diameter). It is for impact conditions below the Bernoulli limit, 
where elastic-plastic response dominates, that material hardness values are rele ant. 
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The above example explains why ceramic materials are useful as armor against jets and LRPs. 
The standard equation (Equation 24) indicates that normally the target material should be of high 
density to defeat an LRP. Yet most ceramics are of low density. It is the unusually high value 
of for ceramics that explains their use, which is due to the high ratio of compressive yield 
strength and low density. Most ordinance is of lower velocity than these critical speeds. 

One technique for understanding the responses is to use what are known as P-u (Pressure-particle 
speed) curves. Such curves are a well-known tool for 1-D impact solutions. The impact stress 
is given by the intersection of the material Hugoniots. For most solids (except porous bodies) 
die release paths are very similar to the Hugoniot, resulting in paths which are essentially "mirror 
images" of the Hugoniots. To understand the responses for "chunky bodies" we introduce the 
concept of Bernoulli adiabats, by analogy to die Hugoniot being also known as the Shock adiabai. 
The Bernoulli adiabat applies only to compressive paths, with release following die normal paths. 

This approach immediately explains an interesting anomaly: for 1-D impact of a low impedance 
projectile into a high impedance target die projectile rebounds after the first stress wave 
reverberation through the projectile, yet for a chunky projectile (same materials and impact speed) 
die projectile does not rebound. The explanation is as follows. 

Immediately upon impact the pseudo 1-D stress is given by the Hugoniot state. When die lateral 
stress relief reaches die impact axis the stress drops to the Bernoulli stress. If both materials 
undergo Bernoulli flow this stress is given by the intersection of die two Bernoulli adiabats. If 
only one material flows then that material jumps to its Bernoulli locus but the other material 
remains on its Hugoniot. If neither material flows then both materials remain on their Hugoniots. 

Figure 1 shows die case of the low impedance projectile impacting the high impedance target. 
The initial stress rises from (A) along die projectile Hugoniot to die intersection with the target 
Hugoniot at (B). This stress then rapidly drops to the Bernoulli state (D). The rear surface 
release wave now originates from the Bernoulli state instead of from die Hugoniot state (path DE 
instead of BC). Thus die projectile velocity merely decreases instead of reversing. Subsequent 
recompressions follow the Bernoulli adiabat, while release follows die normal path. Thus a series 
of reverberations in the projectile occur as the latter comes to rest. Because the projectile is 
sim ultaneously thinnin g as it spreads laterally each subsequent reverberation transit time becomes 
progressively shorter. Thus the projectile rapidly stops. 

The same logic also explains the very deep craters formed when "strong rigid" projectiles (e.g., 
tungsten carbide) impact "weak soft" targets (e.g. aluminum, lead) at relatively low speeds. If 
die impact speed is above u** for the target, but below u^ for die projectile, only the target will 
j um p to its Bernoulli locus while the projectile remains on its Hugoniot, as shown in Figure 2. 
The projectile slows down via a large number of reverberations, and since the projectile does not 
si gnifi cantly deform these reverberations have constant transit time. Thus the total time taken 
to stop is long, and die crater is deep. 
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This d articular logic also applies to the case of trying to capture a particle totally »ntact (no flow 
except elastic-plastic deformation). To do so requires the impact speed to be below the ««* o 
the particle but above the u^ of the target Since the particle properties are fixed the only ophon 
is tod^ose a target with a very low Uerit value. Referring to Equation 3 1 we see mat this implies 
the need for a target material with very low yield strength and vejy low density, ihis des^Ph°n 
applies to the case of trying to capture micrometeoroids with very low density foams and 

aerogels. 

He result of the strength terms is to produce signifies. ^ imp “’ 

speeds. Thus none of POD's equations are rigorous for speeds below the 1 km j range. 

23J& General Comments on Scaling (Similarity) Laws 

The above approaches reveal that it is rarely self-evident as to whether the analysis should be 
JLdt, . momentum logic or tut energy one. His dilemma has important con s equeiK« ance 
momentum conservation logic (as used above) gives an inverse-square .aw (< =•. °, « ■) 
sness amplitude versus propagation distant*. whereas «ergy ^rv«mn ro s gesBa .ample 
inverse law (o oc r*'), for spherical divergence. Reality suggests that the try - response is a 
mixture with m index between -1 and -2. Further, the known non-linear, ties of hydrodynamics 
suggest * that the index is probably variable also, depending on stress levels fe nd propagatoon 
distances. Such effects explain why the various power indexes are difficult to -enve, and 
im ply that these indexes may not even be constants. 

An obvious question is whether the process of impact cratering should obey a 
law Here are reasons to believe the answer is rigorously no aldioug a 
ap proximation " may exist He basis of this comment rests on work done by r ich scienhsts as 
Zeldovich (1992) who studied the following simple case. 

Imagine a semi-infinite half-space of gas of initial finite uniform density hut of aero 

and pressure (this problem is clearly a "thought expenment since such a gas state is 

Allow a piston to impart a transient 1-D push to the free surface of the gas for a Jiort period and 

then be rapidly withdrawn faster than the release motion of the gas surface, th is transferring 

ballistic impulse to the gas. 

He result of this impulse is to propagate a shock into the gas and simultaneous* cam* the free 
surface of diegas to expand outward. Hus a "turning point" exists within the gas dmthng the 
modern ^ions of inward-going mid otnwatd-going, mid dus p, m, ..self follows 

die shock front but at a lower speed. 

Zeldovich attempted to find a single similarity rule that described die behavior Jd, in Partic ula r, 
predicted the rate of decay of the shock pressure versus propagation distance He immecU at y 
, problem' it was impossible to establish a single role that simultaneously conserve ** 
rmenmm mid energy. Cover, Zeldovich was able to find - ^ 

both conservation rules provided a piece of the problem was ignored. Specially, g 
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of gas at the leading edge of the free surface expansion was ignored then the remainder of the 
problem could be described by a scaling law which satisfied die conservation laws. This "slug" 
contained the gas mass which was initially compressed during die ballistic impulse. 

Thus Zeldovich concluded that a scaling law could be asymptotically deter min ed as the "slug" 
mass became negligible compared to the total mass in the moving gas, which itself continuously 
increases with time. For the case of a diatomic gas (y= 7/5, ratio of specific heats) the result 
was that the shock pressure decayed with the 4/3rd power of the distance. Note that Zeldovich's 
problem was simple in the sense that it was 1-D and involved no material strength or phase 
change effects. Note that die power index was between -1 and -2, even though this is a 1-D 
condition. 

POD is of die opinion that a similar situation applies to cratering. It can be anticipated that die 
solution is really one which asymptotes when a portion of the problem is ignored. This portion 
is probably related to the projectile mass and the mass of the target involved in the initial shock 
wave phase. Thus the solution may have two parts, one describing the initial material response 
and one describing the later cavitation phase. As the mass of the cavity (i.e., crater) becomes 
large compared to that of the projectile and initially compressed target material so the response 
can be expected to approach a condition that can be described by a scaling law that satisfies 
conservation rules. This logic implies that simple scaling laws are to be expected only for impact 
conditions which produce large craters (i.e., high speed, low target yield strength, high ratio of 
projectile to target densities). 

2.3.9 Supralinearity 

Although supralinearity has been experimentally observed it is not predicted by any known 
analysis based on hydrodynamics, neither do purely hydrodynamic codes predict it. Suggestions 
have been made that the effect is related to the fact that the projectile flattens (i.e., "pancakes") 
as it impacts the target, and that this modifies the apparent diameter of the impactor. This logic 
has then been coupled with a suggestion that small projectiles are "stronger" than larger ones, and 
thus deform less. 

The hydrodynamic arguments given above involved the term r p 3 , this being a product of the cross- 
sectional area of the projectile and the effective "time-width" of die projectile. In reality, the 
term r p 3 is a measure of the volume (or mass) of the projectile. Since the latter remains constant 
during the initial impact the logic suggests that the exact shape of the projectile is not important, 
at least for small changes. This accords with die common assumption that spheres can be 
approximated by the equivalent-volume (chunky) cylinders and that cratering is a function of the 
projectile mass. If this logic is valid, then the "pancaking" process should have no significant 
effect. Thus supralinearity must be caused by some other mechanism. 

A phenomenon that could explain die supralinearity is die "Petch Law" (Petch, 1953) which 
describes the strength of a ductile material versus its mean grain size. This law states 
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Y, = Y* (1 + (A/d )'*) 


(33) 


where Yq, is an intrinsic yield strength, Y, is the observed strength, d is the me^ii grain size, and 
A is a material-specific "grain size" parameter. The law predicts an increase ; n yield strength 
as the grain size decreases. Experimental data for aluminum (Anderson, 1990/ s ugg ests that A 
is about 50 microns. Petch's theoiy involved consideration of the shearing tha occurs between 
grains. Based on this theory the quantity A is given by 

A = 7tGu/Y a 2 (34) 

where G is the material shear modulus and % is the material surface energy per rea for opening 
cracks. 

POD thus suggests that supralinearity is related to target strength properties, an . that this effect 
is related to crater size and projectile size. However, only the former is important because of the 
arguments given above. Specifically, the important aspect of the problem is tha ratio A/r e , die 
ratio of material-specific size and the crater radius. 

To understand the logic, consider that the normal "bulk" yield strength of a mat-rial is really an 
average over a volume large compared with the dimension A. If die "active volume under 
consideration ( e.g the crater) becomes comparable to this dimension then the normal bulk 
properties no longer apply. Clearly, if the crater involved only a single rmerial grain the 
governing yield strength would be that of the grain rather than of die average material. On the 
other hand, if the impact were on a grain boundary the apparent material strength would be 
lower. Since the probability of hitting a grain is higher than hitting a boundary for a small 
impactor die average response will be a higher yield value. 

To apply die Petch Law we use Equation 33 to define die effective yield strength in Equations 
21 and 28, and replace the term d with r e . Unfortunately, the solution requires iteration of the 
resulting equation owing to the form of the r e terms on either side of the formula However, the 
asymptotic solution can be readily obtained for A » r c . 

Let = (A/YqJ / (1 + (A/r e ) ,/2 ) =o(A/Y<J (35) 

with the term A containing all the details of density, impact speed, etc. Then he asymptotic 
solution becomes 

(r</r p ) (N * 1/2) =(A/Y 0 J(r p /A) ,/J , or (r^jnr/Ar " (36) 

Note how we converted a dependence on r e to one on r p . Thus if die effective jower index is 
N = 3.5 (equivalent to the index used in Equation 21, which employed N+l ratl.jr than N), we 
obtain 

( r e/r p ) = (A JYj* (VA) 1 * (37) 
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giving a supralinear tenn with a power index of 1/6 (0.1667), which is thus the -small-limit- 
index. However, the most-quoted supralinear index is that given by Cour-Palais (1985) with a 
value of 1/18 (0.056). It is interesting to see whether the Petch law can explain the Cour-Palais 
quote, without invoking "absurd" values for the power index N used in Equations 35 and 36. 

The use of the Petch law produces a "downgrading" factor, F, in the crater sizes otherwise 
predicted, so we multiply Equations 21 and 28 by this factor 

i.e., F = 1/(1 + (A/r e ) ,/J ) N ( 38 ) 

Using A * 50 microns and N = 1/3 we compute F versus the "normal" crater prediction, r c . We 
find that for r e = 1.0 micron F = 0.499, for r e = 10 microns F = 0.6761, for r e = 100 microns F 
= 0.8367, for r e = 1 mm F = 0.9349, and for r c = 1 cm F = 0.9775. We immediately see that the 
"effective" supralinear index varies with crater size (hence projectile size). For very small craters 
« 1 micron) the index approaches 1/6, while for large craters (» 1 cm) the index 
approaches zero. Thus the apparent index, n, must be in die range 0<n<l/6. What is relevant 
is that comparing values over a range of crater sizes can lead to an apparent index in die region 
of 0.034 to 0.073. Specifically, over the range 10 microns to 1 cm (or 50 microns to 1 mm) the 
apparent index is 0.0534 (or 0.0546), which is very close to the Cour-Palais value. Thus the 
Cour-Palais index may merely be an "apparent" index which applies over a restricted range of 
projectile sizes. Figure 3 plots F versus r e . 

The corresponding value of A for Teflon is not well known since no quotes have been found for 
the x term. If this term were similar in value to that of aluminum the value of A would be a 
factor of about 1.5 larger than for aluminum taking into account die lower yield strength and 
lower shear modulus of Teflon. This would have die effect of increasing slighdy the 
corresponding values of crater size (rj required to produce a given F value, which in turn means 
that die size range over which the supralineanty is observed would increase slighdy. 

If die above logic using die Petch law is correct then a veiy important corollary results: the 
supralinearity does not apply over all crater sizes but is in reality a small- size decrease in crater 
dimensions. In particular, die effect disappears for craters above about 1 cm diameter. It should 
be noted that most (all?) of the experimental data illustrating the supralinear effect have involved 
projectiles in the size range of microns to mm, which is consistent with the above logic and 
expectations. 

Two other factors could also produce a supralinear effect, namely thermal conduction and strain- 
rale effects. If the strongly shock-heated regions were to rapidly cool via thermal conduction to 
deeper lying regions then the operational yield strength would increase. Normally, thermal 
conduction is ignored in hydrodynamic calculations owing to die fact that stress waves move 
much faster than thermal diffusion waves. However, this difference in propagation rates is small 
for very small scale geometries (below a few microns for metals). Thus very small craters may 
become thermally "quenched" and grow less than would happen without conduction. This effect 
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would be strong only for the smallest craters, and should be negligible for cr ers above 10 
microns. 

,„ lf h « a.-, -reposed as an important mechanism for supralineanty, based on the 
a „ .l I j strongth of many materials noticeably increases for very high rate 

SSs Howevet'titis 2? has », been well proven in hydrocode ans.ysts. Further, 
conditio . , th t t u e craterinc is a direct function of the she k front. The 

increase^n yield stt^gth is observed for the shock front, but since the effect is related 1 to 
flow it is also observed that the stress relaxes back toward the normal yield V5 “ e ^ , e 
S front Thus if cratering is related to the integnd of the str^s pulse versus nme the stram- 
rate effect will only slightly modify the cratering versus the absence ot the effert. A more 
^dfSTproWeL is Ae^ fact that strain-rate h as the unih i of time, -t l^ f T^obtmn a 
lensth narameter would require the combination of velocity/strain-rate. Indeed, if strain-rate is 
responsible for supralineanty then there should be a velocity dependence, such the: higher impart 
S give ZSL "smaller" craters than would normally be expected Tlu. would reved 
S.,f as an ^parent drop in the velocity index. Proponents of stre.n-r.te log.c never mentton 

this fart. 

23.10 Brittle Material Response 

Because TFE Teflon can behave in a brittle fashion it is appropriate to briefly discuss such 
f^^nie major difference between a brittle material response and that of a ductile one is 
that the former can produce radial cracks (centered on the impact point) and an cuter sp cr x 
surrounding the normal yield crater. Brittle materials fail readily under simple tension by te 
“lowever, under compression (especially a tri-axial state with simultaneous latera^ 

toZmtioos), while alumina has a yield strength of over 60 kbars (as evidenced by « dtsplaymg 
a Hugoniot Elastic Limit under 1-D shock compression). 

We can develop an equation describing radial cracks by equating the hoop si ess logic given 
^rr^necess^y to «us« emek gmwfl. in dte mode-1 nmnner, Ibts larter sness 

is given by 


<W = « K i«/ r 


.1 n 


(39) 


where K,. is the fracture of ™ 

“To 
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spaced around the impact point means that the "far-field" distance can be no more than the half- 
distance between these cracks. Thus a is best determined by empirical experiment. 

Equating the hoop tension (Equation 8) with Equation 39 gives 

r cnd/ r P = (P.fl/aK le ) 04 Uo 01 r “ (40) 

where Prir = p p p/(p p 1/2 + p^) 2 (41) 

The following should be noted: 

(1) The equation predicts supralinearity with an index of 0.2. Cour-Palais (1985) quotes 
this same index as being the one for brittle material response for the spall craters versus 
the index of 0.056 for ductile materials. 

(2) The velocity index is 0.8 rather than the roughly 2/3 value for ductile craters. 

(3) Because of these differences in the two indexes it follows that "cross-over" points 
exist. Keeping all other factors constant, if the variable is either velocity or projectile size 
there will be a critical velocity or critical projectile size where r^ = Projectiles 
which are both large and at high speed will tend to produce radial cracks which extend 
out further than die corresponding ductile craters. Materials with low yield strength and 
high values of fracture toughness (e g., most metals) will tend to have r^ < r For 
this case the cracks are "lost" within the ductile crater. For materials with low fracture 
toughness and high yield strength (eg., most ceramics and glasses) the reverse occurs. 
The ductile craters are limited in size while extensive radial cracking will be observed 
extending beyond the crater. 

Experimental impact data from H6rz (1992) on TFE Teflon indicate spall cratering surrounding 
a yield crater. However, when the samples are sectioned (see Figure 4) a series of radial cracks 
are observed. These cracks are roughly equi-spaced around the hemisphere centered on the 
impact point. The hemispherical surface which circumscribes these cracks is observed to intersect 
the original target surface at the radius of the spall crater. Thus the spall crater radius is 
intimately linked to the radial cracks, and the latter define the spall crater radius. The evidence 
implies that the spall crater depth is determined by "peel-back" and fracture of the outer portions 
of the material between the original surface and the first radial crack from the surface. Thus 
Equation 40 gives the spall crater dimension. However, at present there is no simple manner for 
predicting how many radial cracks will be produced, and thus how deep the spall crater may be. 

23.11 Phase Changes and Momentum Enhancement 

Thus far the phenomena of melting and/or vaporization have been discussed only in terms of the 
induced fluid state which has negligible yield strength, shear modulus or tensile strength 
However, as vaporization starts to occur another factor must be included, namely the 
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enhancement of the impact momentum caused by the backstreammg of the vapor Although there 
is always some backstreaming material during cratering (e.g., the jetted material and parts of the 
crater lips which break off) the momentum associated with this is usually small if the material 

states do not include vaporization. 

The major effect associated with vaporization is a dramatic change in the behavi r of the release 
adiabat (isentrope). For solids and liquids the release adiabats are steep and nave very little 
curvature in the pressure-volume (P-V) plane. Thus the release wave speeds arc high and there 
is only a limited spread in the wave speeds (wave speeds are given by c - -V dP/dv, where V 
is specific volume). Furthermore, the expansion required to reach zero pressu e is also small. 
Consequently, the stress relief is "rapid". The released energy (i.e., PdV) conve ed into motion 
is also small (i.e., the shocked material remains hot even at zero stress). 

For the vapor state, however, the release adiabats have significant curvature in th P-V plane and 
significant expansion is required to reach zero pressure. Hence the average wav speed * lower 
and there is a significant spread in release wave speeds giving the well-known "Taylor Fan The 
time to reduce die stress is much longer than for the solid/liquid states. Thus the extra 
momentum is a consequence of the larger release pressure-time integral for vapors. The released 
energy for a vapor is large and ultimately becomes converted totally into kinetic e tergy of motion 
and thus the temperature adiabatically decreases as expansion occurs, for a full vapor state (i.e., 
one whose adiabat passes above the material Critical Point). 

For release states which cross die two-phase liquid-vapor condition there is £ gradual change 
from die liquid behavior to the pure vapor behavior as the initial shock loading increases. 

If cratering is dependent on net momentum then the vaporization will cause deeper craters, but 
will have little effect on crater diameter since the latter is mosdy determined uoring the initial 
reverberation within the projectile during which the "long release tail" has not developed^ For 
total vaporization of the projectile the net momentum increases by considerably less than a factor 
of two, partly because the vapor state of the projectile occurs upon generation of the first release 
wave from the back of the impactor, and at this time the projectile is moving fo; ward with about 
half the initial impact speed, but also because the vapor cloud expands hemisphencally upstream. 

Since half the projectile momentum (for a like-on-like impact) is transferred to Oie target during 
die first shock wave motion from the projectile/target interface back to the target free surface, 
only the second half of the momentum can "bounce". Thus at worst we expect the total 
momentum to be 1.5 the initial projectile value. This is an upper limit since it assumes the 
rebounding material remains as a solid particle with no lateral spreading,. 

An approximate evaluation of the additional momentum for total vaporization can be obtained 
by using the logic developed by Gurney (1943). The essence of Gurney's theories is that die 
expanding vapor cloud rapidly asymptotes to a condition where there is a leading edge velocity, 
„ which becomes a constant, and the gas between the source of the gas expansion and the 
leading edge, at r. , has a linear distance relationship of the form u(r) - ru„ t m . The second 
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assumption is that, at any time, the density in the vapor cloud is a constant determined solely by 
the volume swept out by die advancing edge r„ versus die source. Using these rules, and 
assuming hemispherical expansion from a fixed source we integrate the energy wi thin die 
hemisphere and set it to one-half of die impactor energy (since one-half is transferred during die 
same wave motion) to get 

u.« (5/6) 1 * uo (42) 

We now integrate over die hemisphere to obtain die net directed momentum (i.e., along die 
impact axis) to get 

I = (3/16) (5/6)** Io (43) 

where Iq is the initial projectile momentum. Thus die total momentum becomes 

Imi = 1.171 ^ (44) 

A further refinement occurs when we assume that die vapor expansion occurs from a mov in g 
source. For a like-on-like impact this speed is -u</2, where we define a positive speed as the 
backstreaming motion, and the initial impact speed as negative. Hence the velocity distribution 
becomes 


u r = -u 0 /2 + r/r„(u„+u 0 /2) (45) 

Note that this implies there is a "turning point" within die vapor, which divides the flow into two 
regions, one following the projectile/target interface and die other flowing back upstream. Next, 
for a like-on-like impact one-half of the kinetic energy is transferred to die target during die first 
shock wave transit through the projectile, before die vapor expansion starts. We therefore equate 
the remaining half-energy with the energy in the hemispherical cloud by integrating 


0.5 KE p = 7tr p J p p Uo 2 /3 » J 2%tp t (u^/2) dr 

(46) 

where p p = 2p p (r/r. ) 3 

(47) 

Solving, we obtain: u»= +u</4 

(48) 


Thus the vapor backstreams with a leading edge speed of one quarter of die initial impactor 
speed. Integrating over the hemisphere for the momentum we find 

Momentum = I</32 (49) 

Thus die total momentum becomes 


= Io + V32 = 1.03 Io 


(50) 
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Thus the effect of strong vaporization of the projectile is to produce an increase in the total 
momentum bv a factor in the range of 3% to 17%, with the true value closer to tl : smaller quote 
STTlirt-to generates much of its motion before tit. projectile amps, for tins reason 
the effect is not readily identified in terms of crater sizes, and die cratering response will appear 
to be "smooth" with no "jumps" as vaporization occurs. This implies that ulti -hypervelocity 
impact responses will appear to be continuations of the lower speed impact responses. 

Because a portion of the target crater can also be vaporized the momentum « nam* 
more However, since the peak stresses drop rapidly with propagation distance th target material 
subject to vaporization is usually only a small fraction of the crater ; exce P t for very high imp* 
speeds. However, it is necessary to be aware of a "chicken and egg probiem. Lawmw (1989) 
has predicted a significant momentum enhancement for very high speed impacts The soiirceof 
this enhancement is the large mass ejected from the crater which is assumed to m m the ^vaj>o 
state and to contain a significant fraction of the projectile kinetic energy However, is the mass 
ejected greater because the momentum was higher, or is the momentum higher because th 
ejected is greater? POD believes the latter is the correct interpretation. If so, the momentum 

. larger net push to tire urge, bu, may have Utile effect - ««ermg 
dynamics. This is because a large fraction of this momentum constats of a long-lived low 
pressure state, and this low pressure will tend not to cause much further cratenng. Thus if the 
pressure decreases below that necessary to ensure die Bernoulli state in the target, the cr *'' nn 8 
efficiency rapidly drops. Overall, the effect of strong vaporization of the target is probably to 
deepen the crater somewtwtt, but not necessarily by a significant amount This natter deserves 

further study. 

23.12 Dimensional Analysis 

Many investigators (Holsapple, 1987; Herrmann, 1986, and others) have used the techniques of 
■Dimensional Analysis" to derive cratering scaling laws. Such techniques empioy Buc gb“n 
"pi* groups of parameters, as described in Bluman (1981). Thus because toe quca«.ts (<W and 
(P/ilT are* both dimensionless, toe "pi" groups must likewise be drmensionl«s How *? r ' 
to^is really a mathematics-based one rather than . phjuics-based one 
necessary to correctly identify the material parameters which make up toe pi groups, es 
erroneous answers can ensue. 

For example for densities it is common for the ratio (p/ft) to be used. A con ssponding group 
Led by POD, however, is ((p/Pj/fl-KP/Pj'”)’) which is wuhmPOD’s equations for diameter. 
Whereas the former ratio is "intuitive", the grouping used by POD is not. 

Similarly, it is common to find die groups (pc J /Y) and (u/c) in many seal ig 

however! puts these two groups together to obtain ((p/Y)u ) which is still dimei ^lonless but does 

not include material sound speed, which is handled independently. 

The primary reason for these differing groups in POD’s analyses is that, unlike the Dimensional 
Analysis approach, POD's analyses are based on a consideration of the physic , mvol 
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POD's equations obey the rules of dimensionality as a by-product of the analyses, rather than as 
die starting point for die analyses. POD believes that deriving results on die basis of physical 
logic is a much more appropriate approach than die use of Dimensional Analysis. 

2 A Perforation Behavior 

The above discussion applies to infinitely thick targets. As die target is progressively thinned, 
however, the following sequence of events occurs, assuming all other factors remain constant. 

First, below a certain thickness die target rear surface will display a permanent outward bulge. 
This is due to a region of target beyond die crater depth being subjected to a permanent residual 
compressive plastic stress which attempts to relieve its stress by outward motion. This bulge 
progressively increases as the target is thinned. For slightly thinner targets an internal spall 
occurs (a void) allowing die back surface bulge to protrude even more. However, the bulge does 
not break off, and die void can only be observed by sectioning the target. 

Second, with further target thinning, the bulge splits and lips are formed around a spall crater on 
die target rear surface. This occurs because the reflected compressive stress returns into the 
target as a tensile rarefaction. Beyond a certain distance away from die rear surface, the 
algebraic sum of this tension and the remnant of die still-forward moving compression add up 
to a tensile value which exceeds the local target material strength. This explains the void 
described above. The spalled region detaches itself due to die stored momentum which allows 
die resultant edge stresses to exceed the tensile strength. 

Third, with further target thinning, the final spall surface (there can be multiple sequential spalls) 
approaches the original front-surface crater depth. The resulting relief of stress at the deepest 
regions of die crater also allows continued increases in die crater depth. Thus at a critical target 
thickness the rear surface spall meets the crater depth, and a hole is created through the target. 
This condition is generally referred to as the "ballistic limit". 

Fourth, for even thinner targets, die central hole rapidly increases in diameter. Simultaneously, 
relief waves from the target back surface arrive at die front surface and begin to modify die crater 
mouth, ultimately reducing its diameter. Until this happens, the crater mouth is essentially 
"ignorant” of the target rear surface behavior, and the perforation response resembles a case of 
a "truncated” crater, where only the deepest sections of the crater are affected. 

Fifth, even thinner targets cause die central perforation hole to pass through a maximum in 
diameter, and to then decrease. Simultaneously, the crater mouth diameter progressively 
decreases, as does the rear surface spall hole. These two latter diameters are always larger than 
die central hole. 

Sixth, for very thin targets, all three diameters decrease in a systematic manner, becoming closer 
in value to each other and also to die projectile diameter. Ultimately, a hole is formed which 
asymptotes to the same size (and even cross-sectional shape) as the projectile. Thus for very thin 
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targets (i.e., T<dp/100) the projectile essentially punches its own cross-sectional image through 
the target. This effect is of major use for identifying projectile dimensions, since the target gives 
a direct measurement of the latter, while the damage done to die projectile minim al. This 
minimal projectile damage is a consequence of the very short-lived shock pulse which enters the 
projectile. Such thin pulses undergo very rapid hydrodynamic attenuation. Thu# the stress levels 
rapidly drop and any damage in die projectile is limited to a thin skin on the npact surface. 

All these events have been well characterized by laboratory impact experiment However, such 
experiments have been done over a limited range of impact speeds, or with only very small 
projectiles (eg., micron sized) and with a limited range of impacting projectile and target 
materials. It therefore remains an issue of establishing die physics of the processes in order to 
anticipate die responses for generic conditions. 

2A.1 Stress Wave Logic: The Ballistic Limit 

We shall attempt here to describe the processes leading to the Ballistic Limit ondition. Note 
that the latter term has more than one definition. For some researchers it represents the condition 
of just producing a through-hole; for others it represents the case where the projectile just passes 
through the target (i.e., d,, = dp). However, die change in target thickness (all other factors 
remaining constant) between these two conditions is quite small. Since POD is attempting to 
establish physical models for development of scaling laws we believe that such "nuances" are 
beyond the capabilities of simplified analysis. Consequendy, we shall describe a generic Ballistic 
Limit condition that does not attempt to distinguish these subdeties. 

The logic invoked to explain front surface cratering implies that a spher ;ally diverging 
compressive stress wave moves through die target as if centered on die iir act point. At 
distances beyond the crater depth the momentum law is expected to apply. Th ; gives 

= o 0 (rp/r) N (51) 

where q> is a source stress originating at the projectile/target interface, r p is the projectile radius, 
r is the local distance from die impact point, and the index N=2 for a non-dispers /e system. We 
assume that the source stress is the Bernoulli one, given by 

q, = 0.5p p u o J /(l+(p p /p k ) ,/2 ) 2 (52) 

and we also assume that when this stress wave reflects from die target rear surface it continues 
to spherically diverge as if centered on an image point at twice die target thickness on die impact 
axis. Thus the stress amplitude inverts into tension, but otherwise acts as if the wave continued 
to move into matter beyond the back of the target. The algebraic sum of this tensile wave and 
the remains of the forward-going compressive wave gives the net stress. The forward-going 
compressive wave drops its radial stress to zero at the crater surface. Thus the reflected wave 
generates the maximum tension as it approaches the crater surface. 
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Assume that if this tension exceeds the local target spall (fracture) strength, o., die material 
fractures. If this fracture occurs close to the crater bottom, then the total propagation distance 

for die pulse was 

r = T + (T-P) = 2T-P < 53 ) 

where T is die target thickness and P is the crater depth. Thus we rewrite 

((2T-Pyr/ - Cfe/o. , or (2T-P)/r p = (oM” < 54 > 

But (2T-P)/r p « ((2T-P)/P)(P/r p ) ( 55 ) 

Hence 2T/P = 1 + (r p /P)(q/o.) ,/N (56) 

and T/P = (1/2) (1 + (d p /2P)(<yo.) ,/N ) (57) 

Alternatively, since T/P - (T/dpXdp/P) we may rewrite to obtain 

T/dp = 0.5(P/dp) + 0.25(q/a,) ,/N ( 58 ) 

This equation thus states that there are two components to the condition causing the spall plane 
to coincide with the crater bottom. 

Substituting Equations 28 and 52 we obtain 

T/dp « (l/8)(4/3) l/3 (p p /p I ) 1>3 (p l /YJ 1/S { (c* t +s(u,-ii Vra *)/( l+(p t /pV) 1/I ))(u,-u Vc J} 1/3 

+ (l/4){ p p u # 1 /(2ai(l+(p l /p l ) l/2 ) 2 )} 1/N 

If N=2 the second term becomes 

(1/4) (p J /2o.)' /2 Uo /(l-Kpp/pk) ,n ) (60) 

Several relevant points can be made: 


(1) The equation contains two independent parts. 

(2) The equation employs both yield strength (with a 1/3 index) and tensile strength (with 
a 1/2 index, if N=2). We believe this explains the confusion over which strength term 

to use for perforations. 

(3) The two parts of the equation have different velocity indexes. The first term has an 
effective index that starts below 2/3rd but asymptotes to 2/3rd at high speed. The second 
term has an index of about 1.0 (if N=2). Thus the combination of terms will appear to 
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have . velocity index (M) in die range 2/3 < M < 1.0. 

itself not be consunt, but will slowly charge Ae BJMc 

decreasing if N>2. All previous existing equations (see Secti ) 

Limit have a velocity index of 2/3 < M < 1.0. 
o , . nw _ the effects of variations in N on Equation 59 for aluminum on aluminum impacts. 

approximation for aluminum. 

T , . . . nflted u.-* as the Ballistic Limit condition is approached, the reflect stress wave 

It should be noted that as me mi development If the crater 

Z^gTholt b How'e?er! LZc«* crator has no, completely sjopped M spall occurs, 
Ae continued hoop stretching causes this thin wall to open up into the hole. 

The CTH data indicated that the deceleration of die crater bottom was at a pseui -constant rate. 
£Le™ e meTs^eed was almost one-half of die UH speed genemtod a, im,act 

Thus u„ = 0.5 u/d'Kp/P,) 1 '’) ^ ’ 

and the time to form the crater depth is 
x* P/u„ 

During this time the stress waves moves at speed c, and propagates a distance 
X = C t P/u„ 

Thus x/P * 2 c, 


(62) 

(63) 

(64) 


For an Al/Al 1 impa«, , dtisratic > ^ Slow 

cttor before die latier has iully developed. 


v/P = 2 (Crw + S Uo/d-KP/Ppi'^Kl-KP/Pp) 1 ^ 




which for high speed impact reduces to 
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x/P * 2 s 


( 66 ) 


For aluminum s= 1 . 4 , therefore x - 2.8 P. But the required propagation distance for the stress 
wave to reflect and reach the crater bottom is x * 2T - P. Thus, T = 1.9 P is the maximum 
limiting target thickness (for aluminum) at high impact speed that allows the reflected wave to 
pass the crater bottom before the latter has finished its full "normal" (j.e., infinite target) 
development. The above quote should be compared to die quote of T = 1.8P to just prevent 
perforation, as given by Christiansen (1992b). 

For TFE Teflon s = 0.795, therefore x = 1.59 P. But the required propagation distance for the 
stress wave to reflect and reach the crater bottom is x = 2T - P. Thus, T - 1.295 P is die 
maximum limiting target thickness (for TFE Teflon) at high impact speed that allows the 
reflected wave to pass the crater bottom before the latter has finished its full "normal" (i.e. 
infinite target) development. 

2.4.2 Ultra-Thin Targets (Foils) 

For perforations of ultra-thin targets (^d/lOO), we assume a similar logic to that used to 
describe crater diameters. Again, we assume that the radially propagating pulse contains 
momentum which induces hoop strains and stresses, and that the cratering process stops when 
die latter drop to the local yield strength of the target. Thus 

Y» = ajdt c/r e (67) 

However, there are two major modifications for ultra-thin targets. First, the effect of rapidly 
arriving release waves from both the front and rear target free surfaces causes the stress to 
decrease with propagation distance very rapidly, such that 

<* = <*> (V r ) N (68) 

with N » 2 . Second, die effective time- width of die pulse is no longer related to the 
reverberation time within the projectile, but instead is given by die reverberation time across the 
target, since beyond this time the projectile has effectively punched through the target and there 
is no further lateral push. Thus 

dt = 2T/c, ( 69 ) 

Thus from Equations 52, 67, 68 and 69 we obtain 

Y» = (Pp/2)( 1/(1 +(Pp/ P t) 1/J ) 2 Uo J (2T/cJ (c/O (r,/r e ) N (70) 


which gives 


( r </ r P ) N = (Pp Uo 2 /O+fPp/ptD 1 } (T/r e )/Y, 


(71) 
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» (72) 

the values of 
nt point, since 
asymptotes to 
'een the crater 


of it becomes 
drive" on die 
ion due to the 


w<v - wg - (p/py^^P« lw *'x w ‘ ,, P T ^ w * , v(<^prr 

Pl l j» m » 2 this equation collapses toward ^ 1 -P’ tegasdlessof 

S^H« v,eld stren«h target thickness or impact speed. Hus latter nt an .rnporta 

JS-SSSS-- observance. Thus the “““ “ fferer' “lS 
the same size as the projectile. For this case there is no significant difference 

"mouth" diameter and the central hole diameter. 

_ . . r , xt •» 9 ;« that although local momentum must be conserved, muci 

e,C H £ -d££ and back lips and, therefore, the -effechv 

«*-*• Al “- ”*" ow pulses und “ go rap,d 
release waves which overtake the shock front. 

2.43 The Intermediate Case 

Most trackmg 

ssr «\7i - 


N = 2 (1 + m d/T) 


(73) 


where m is a multiplier whose value is asto mdce onlhe^Ttha? the 

cratering equation, Equation 21, when T -2/3 con ^tion. Thus N 

crater diameter is almost identical to the in * ni ® . to 2mdJT for T « d. for ultra-thin 
asymptotes to 2 when T » d, (the mfirnte target ca «), and to 2n«vr u, 

foils. P Further, dividing Equation 72 by Equatton 21 we obtain 

(74) 

(d^cipf- 25 - 3T/2dp 

which can equate if N = 2.5 when T = 2/3 dp. 

Thus we chose m * 1/6 - 0.16667. Hence 

N = 2 (1 + 0.1667 dp/T) 

... , /nu+n = o 2857 for Equation 72, which is the same index used 
?£££«• ^StISb £• use Equation 2., white hrKM d,w. u* Equmion 72. 


(75) 
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Figure 6 plots Equation 72, using the index in Equation 75, for (d/dj versus (<L/T). Also shown 
on this plot are the experimental data from Hftrz for (<ydp). We see that tS trends are very 
amilar POD does not claim that the final equation has been derived; rather we believe that the 
fonn of the behavior has been identified. 

Allying the logic of Equation 75 within Equation 72 also leads to the following behavior as 


(d, - dpVdp = (3.0 T/dp) ( In (A) - In (T/dp)) 
where A - 2 (Pp/ Pt ) (p/YJ ^ / (1 + (p p /p l ) ,/3 ) 2 


(76) 

(77) 

This suggests that the hole closes down to the projectile size in a roughly linear fashion as taraet 
,s reduced for vet, due ttrgets. However, the loganthmic term for (T/<L) pr^2“ 
net effective index which is somewhat less than unity, especially as T/(L increases For 
example, for 0.005<T/dp<0.01 the apparent index is about 0.894, while for 1.0<T/<L<2 0 the 
apparent index is about 0.776. These index quotes should be compared to values ofabout 2/3 

K*.* S *Tl e „l I96 2: K ? id “ (1963 > md Brow " O 97 ®). •”<! die quote of 0.895 given by 
SAonberg (1988). The loganthmic behavior for the term A demoestra.es a ve ' weak 
dependence on material properties and impact speed, in accordance with observations. 

e f e f of variations in N on Equation 59 for aluminum on TFE Teflon 
A * f j- . calculations done for TFE Teflon cratering in infinite targets demonstrated 

that, for distances greater than the crater depth, the stress decay roughly obeys a 1/r law near the 
impact point but steadily asymptotes towards a 1/r 24 law with increasing distance Thus N=2 4 
is a good approximation for TFE Teflon. LA 

Starting with Equation 72 we assume that the index N is itself a function of T/d As for 
aluminum a simple possibility is the function ^ 


N = 2.4 (1 + m d/T) 


(78) 


^, Ch . h ^” 9 ' 0t “ t0 2 4 . wh “ T » J (die infinite target case), and to 2.4md/T for T « <1 for 
uhra-Am fods. The value of m is chosen so as to make Equation 72 identical to Equation 21 
when T - 2/3 dp. Thus we chose m = 1/36 = 0.0278. Hence 

N = 2.4 (1 + 0.0278 d/T) (79) 

Apptymg the logic of Equation 78 within Equation 72 also leads to the following behavior as 


where 


( d « * dpVdp = (15.0 T/dp) ( In (A) - ln(T/dp)) 
A = 2 (Pp/ Pl ) (p/YJ Uo J / (1 + (Pp/ P J ,/2 ) J 


(80) 

(81) 
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2.4.4 Perforation Supralinearity 


Experimentally observed supralinearity for aluminum targets was explained oy invoking die 
Petch Law" (Petch, 1953) which describes the strength of a ductile material versus its mean 
grain size. POD anticipates the size effect to also apply to perforations, and to tause die critical 
target thickness for ballistic perforation to scale with projectile diameter in a similar manner as 
does the crater diameter. Thus the above equations describing die Ballistic L.mit, etc., should 
be "downgraded" for very small projectiles, by a factor F, as described in Section 2.3.9. 

2JS Oblique Impacts: Summary of typical Behavior 

All of the above discussions referred to impacts normal to die target surface. Tie following is 
a brief description of typical responses for oblique impacts. 


Defining the angle 0 as that between the projectile motion and the normal to ti ; target surface, 
it has been observed in experiments that for 0 < 60 degrees the major effect for the craters 
to develop as if the impact speed were given by i^cos©. Therefore, all of die existing scaling 
laws are adjusted for oblique impacts by adding in a cos0 correction. Thus as the component 
of the impact speed normal to the surface decreases so the crater also decreases -i both diameter 
and depth. Further, the crater aspect ratio changes very litde, and die craters art essentially still 
axially symmetric. 

However, it has also been observed that for 0 > 60 degrees, the craters start > become more 
asymmetric. The "downstream" portion tends to become elongated. At this stage portions of the 
projectile also start to ricochet and material is ejected downstream at a small angle to the target 
surface (typically within about 15 degrees). Other phenomena include observations of the 
projectile itself shearing, such that the upper portion (away from die impact surface) can detach 
itself and impact the target downstream as a separate impactor(s). 

2.5.1 Idealized Theoiy 

The general response for oblique impacts can be understood if we start with an "idealized" 
condition; namely, the assumption that die interface between the projectile and die target has a 
zero coefficient of friction. We also ignore the fact that the projectile peneti ,.tes (i.e.,"digs" 
itself) into the target. Under these circumstances, an oblique impact will transmit to the target 
only the component of momentum and energy that is along the target surface normal. The 
remaining momentum and energy reside within the projectile as components parallel to the target 
surface (we call these latter components the orthogonal ones). Since hypervelocity impacts 
produce an effective coefficient of restitution of zero (i.e., there is negligible bounce) there are 
also no components remaining along the target surface normal. Thus the projectile ricochets by 
skimming along the target surface in the downstream direction. 
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The above logic implies that the target response will always be that for the normal component 
of the impact speed only, and also implies that the projectile always ricochets across the target 

surface. 


Reality of course, is a little different. Because the projectile impact immediately causes jetting 
and the onset of lip formation outboard of the projectile diameter, the ability of the projectile to 
simply "skim" across die target surface is suppressed. Consequently, die orthogonal components 
of momentum and energy parallel to the target surface become dissipated by providing a degree 
of asymmetry in the crater diameter. This asymmetry is quite small until the obliquity increases 
to about 60 degrees. The influence of the orthogonal momentum and energy on the crater depth 
is essentially zero. 


Following this logic, POD's cratering Equations 21 and 28 become, respectively. 


P/d,-(l/4K4/3) l “(p;pJ UI (p/YJ“ , l(c.+j(ii,cos9-u u J/(l+(p/P p )‘ , ’)K«,c.^-u^)l" 5 

whereas POD’s perforation Equations 59, 72 and 77 become, respectively, 

T/d, * (l/8)(4/3) lo (Pp/pJ ia (p/Tt) 1 °{(c»t+s(u«cos0-u tia J/(l+(p f /p,) I/2 ))(u t cos0-ii tiO d} 

+ (l/4){pp(u,cos0 ^/(lo.fl-Kp^ 1 *) 1 )} “* 


(d c /d p )*(r e /r p )*(p p /p 4 )* <v,+,) (p/T l ) 1/<N+1> (u,cose )^ +1 >(2T/dp) M<J,+1> /(l+(P J p ( ) 1/1 ) 1 * N+1 > (85) 

A- 2 (jVp.) (p/Y.) (U.COS0) 2 / (1 + (p./ft) 1 ") 1 (86) 

Figure 8 shows a normalized plot of how Equations 82, 83, and 84 vary with the angle of 
incidence (0). The equations are normalized by Equations 21, 28, and 59, respectively. Note 
that the angle of incidence must exceed 30 degrees to get more than a 10% change in T/dp, and 
must exceed 35 degrees to get more than a 10% change in d^/dp or P/dp. In fact, to get more than 
a 20% change (the limit of many experimental measurements), die angle of incidence must 
exceed 40 degrees and 48 degrees, respectively. Also, note that, due to the two components o 
Equation 59, for angles above 80 degrees, T/dp varies linearly. 


2.5.2 The Ricochet Case 

We define ricochet as being the case where die projectile, in whole or in part, escapes from the 
impact crater before the crater lips are formed. For a hypervelocity impact, where both projectile 
and target behave as fluids, a simple condition for guaranteeing that the whole of the projectile 
can ricochet across the target surface is that the orthogonal component of die projectile velocity 
be greater than the speed of the corresponding target surface disturbance. The latter is the sound 
speed in the target, adjusted for the high pressure impact state. Thus we require 
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(87) 


UoSinG^c,,, + s UoCos0/(l+(p/Pp) ,/J ) 


«h.,. s is the Hugoniot term relating shock speed to particle speed, me c„ is he target low 
stress bulk sound speed. Rewriting 

sinG^Co/iio + s cosG/(l+(p/Pp) ,/2 ) 

ITiis equation presets that the entical Juminum into 

, .... *_ -r-neiHiar is the limiting condition which ensures that .he top of the 

lip 'formation 8 This is the limiting erne for p ml projectile 

riwch^t. This castf is defined purely by the geometry of the impact such that 

(89) 


tanG= r/dp - 1/2 (d./dp) 

IT- o rsiratc it«iG versus <L/dL at various velocities. The points where the 2ta iG line crosses 
"T* j /j i; ne i s the limiting angle for partial projectile ricochet for that velocity. The 
corre^ndingcondition for the mid-point of the projectile to clear the crater is g, /en by plotnng 
tanO versus d^dp instead. 

A fund condition to consider is that which ensures the materials behave as fluids. Thus we 
require 

(90) 


UqCOsG ^11^ 


where u„, (defined before) is the minimum speed necessaiy to mam Bemoull flow behavior 
For impact speeds below this value the response becomes elast.c-pl.sbc 

condition. 

r T 7 m «: anp oo oq oo This figure also sh ~ws the two loci 
Fiffure 10 shows a composite of Equations 88, 89 and t> • r a \ir e 

K 1 mid-poin, o?the projectile 

see that for impact speeds in the range o ^ loeic It's interesting to note 

SKCS«=a"R5U — . - 

the Equation 89 criteria used) as is normally seen m experiments. 
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Also plotted is the incipient vapor locus for Al/Al ( i.e ., UoCosG = 10.2 km/s). We see that only 
for veiy high impact speeds above about 30 km/s does this locus catch the ricochet loci. Above 
these speeds the ricochet is of an expanding vapor cloud. 

Christiansen (1992c) has also studied ricocheting effects for Al/AI. For thick targets he derived 
a locus of critical angle versus impact speed. The logic used was that if the stress wave in the 
projectile travelled from the impact point to the projectile rear surface and then back to the mid- 
point before the mid-point itself reached the target surface, then the upper portion of die projectile 
can suffer breakup and cause downstream multiple cratering. Thus defining the mean projectile 
wave speed as c„ (shock out, rarefaction back), the projectile radius as r p and the normal 
component of impact speed as u 0 cos6(with 0 the obliquity) we have 

3 V c p ^r/uoCosG or cosSSc/Siio ( 91 ) 

This locus is also shown in Figure 10. We see that it is close to the loci given by POD. 
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3.0 CTO HYDRODYNAMIC CALCULATIONS 


3.1 Cratering in Aluminum 

The CTH code from Sandia National Laboratory, Albuquerque (SNL4), ha been used to 
investigate the cratering responses of an aluminum target. Parameters varied have included 
impact velocity, projectile density, and target and projectile yield strengths In all cases the target 
was sufficiently thick to behave as an infinite body. For all cases the projectile vas a sphere of 
d«anw»fr»r 100 microns. The resulting data have been mapped and compared to quations of die 
form 

P/dp « constant (Pp/p,) A (p/Y ,) B u, c (92) 

and d^dp = constant tp^ft)/(l +(p/ft)“)*) 4 (fVY^ B u, c (c/Cp) D (93) 

Figures 11 to 20 show die data in log-log form. It is immediately clear tiiat the indexes are not 
constants and that they vary with the impact speed. This effect is strongest at t e lowest speed 
of 1 km/s (which normally would not be considered as hypervelocity) and slowly asymptotes to 
steadier values at the higher impact speeds. 

To date, the best overall fits are given by: 

for P/dpt constant = 0.0725, A ■= 0.60, B = 0.263, C = 0.664 

(Note tiiat these indexes are similar to those in Equation 28, except for the densit index, and the 
leading constant). 

for d^/dp! constant * 0.468, A » 033, B = 0.258, C = 0375, D *= indeterminate (not tested) 

These fits are primarily for data above 5 km/s impact speeds, since die low speed impacts clearly 
are near the limit of applicability of Bernoulli logic and expectations. It is gratifying to see that 
die responses are similar to the expectations. 

Of particular interest are die plots of (P/dp) and (d,/dp) versus projectile density and speed 
(Figures 15 and 16). Note how for die crater diameters the plots "curve over" as projectile 
density increases. This is the form predicted by Equations 14 and 21. For the penetration depths 
die plots follow a more direct power law, which is as predicted by Equation 28. Figures 19 and 
20 show the ratios P/d,. From Equations 21 and 28 we obtain a ratio which Is complex and 
contains low power indexes for density, yield strength and impact speed. For vary high impact 
speed, this asymptotic ratio gives the following equation 

P/d e =>0.253 ( Pp / Pt ) 0<M76 (p/Y J 00476 (c p /c l )°' 2,57 s ,/3 u o 0 w5 ( 1 +(Pp/ft) ,/2 ) 0 ' 57,4 ^(l "KPr P P ) ,/2 ) (94) 
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Specifically, the ratio P/d,, then becomes a positive function of (ty/pt) and an inverse function of 
Y t . The plots show a similar behavior. We thus conclude that much of the basics of the scaling 
laws for either crater diameter or crater depth are indeed contained in Equations 21 and 28, 
although we cannot yet vouch for the exact values of the power indexes. 

Figures 21 to 26 show some of die geometries of the crater for die 10 km/s impact velocity case. 
The dark spots on the impact axis are Lagrangian points, as described below. These figures are 
split, with die right side showing geometry only and die left side giving stress contours as per 
the bar legend to die right of each figure. 

Scrutiny of die CTH runs reveals many interesting facts. It is observed that die maximum depth 
of die crater occurs (in Figure 24) before die maximum crater diameter has formed (in Figure 26), 
and that the ratio depth/diameter is not 0.5 in general (i.e., die crater is generally not 
hemispherical). Various tracer points (Lagrangian points) were placed along the impact axis. 
For these points the time histories of pressure, velocity and position were obtained. Figures 27 
to 32 show such data for the case of a 100 micron A1 ball hitting an A1 target at 10 km/s. Both 
materials were A1 6061-T6 with a yield strength of 2.7 kbars. The Lagrangian points were LI 
through L6 at initial depths of 0.0 (surface), 120 microns, 240 microns, 360 microns, 480 microns 
and 600 microns, respectively. Of these, the points LI to L3 were within die finally formed 
crater which developed a maximum depth of 300 microns, giving P/dp = 3.0. The final crater 
diameter was 500 microns. Thus d^dp = 5.0, and P/d e = 0.6. In all cases the definitions of depth 
and diameter are relative to the initial target surface. 

The fact that these three L points remain on the axis and form part of the crater itself 
demonstrates that the material on die impact axis does not flow out of the crater but merely 
moves into the target. Only target material off the axis takes part in the circumferential rotational 
flow, with that nearest die free surface affected the most. 

Immediately upon impact the peak stress is about 900 kbars. This is lower than the expected 1-D 
stress of 1.63 Mbars (perhaps owing to the artificial viscosity in die code?) but is considerably 
higher than the corresponding Bernoulli stress of 337.5 kbars. This peak stress rapidly drops 
through the Bernoulli value to a very low stress of die order of the plastic yield value, whereupon 
die stress drops much more slowly. Some residual stress reverberation ringing can be seen on 
die plots. The initial shock wave speed is 12 km/s at die 1-D stress, while the release wave 
speeds are even faster. Thus the 1-D stress lasts for less than 10 ns and immediately collapses 
toward the Bernoulli stress. The peak stress then propagates at a wave speed of near 6 km/s, 
which is close to the dilatational speed for aluminum, and decays in value with distance. This 
decay roughly obeys a 1/r law near the impact point but steadily asymptotes towards a 1/r 2 law 
with increasing distance. Thus at early time the stress is obeying the "energy" logic, while at 
later time it obeys the "momentum" logic. The previous suggestion of mixed non-steady power 
indexes indeed applies. 

The induced compressive strain during the 1-D condition is given by (u/W), where u is the 
induced particle speed change and W is the shock wave speed. Thus for a speed u of 5 km/s and 
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a W of 12 lcm/s the strain is 0.417 and the density has increased by the ratio W W-u) vdnch is 
\ T ta Ae concept of treabng the materials as "incompressible’ beta* « com** *>• 
state. Thus the conventional Bernoulli logic applies only after the shock moves away from the 

impact site. 

After the stress peak passes a given point a stress gradient exists between this pulse and the crater 
surface. This gradient decelerates the surface and eventually brings it to rest. Tis deceleration 

is given by 


du/dt = - (1/p) dP/dx 


( 95 ) 


where dP/dx is the stress gradient. If this gradient stayed constant the crater surface motion 
would become a simple quadratic function versus time. However, the gradient anes with time 
and position giving a more complex solution. 

The plots indicate a very small elastic recovery. The stress gradient appears t be attached to 
the rear of the shock pulse at a value of about the target yield strength. 

The CTH calculations demonstrate that the yield strength of the target is indeed a major player 
in the crater formation. Further, a large region (extending to beyond twice the .rater depth) of 
the target is put into permanent plastic strain. This effect is most clearly seen for the deepe* 
Lagrangian points (Figures 31 and 32) where it is seen that the material is ultimately left with 
.permanent compression of about the yield stress. This is the response expected for a material 
which has been subjected to a pseudo 1-D stress wave which exceeds the Hugon it Elastic Limit 
fHEL) Note that HEL = Y(l-u)/(l-2u) where o is the Poisson ratio. For aluminum u 0.33, 
thus HEL = 1.97Y (i.e., 5.32 kbars for A1 6061-T6). Upon shocking the material returns to a 
state of zero stress in the direction of the stress wave but has lateral compressive stresses equal 
to the yield stress. Thus the pressure becomes 2Y/3. At these large radii the ho p stress is very 
small and provides only a small correction. 

A considerable amount of strain energy is stored within this plastic region beyond the crater and 
energy also exists in the propagating elastic wave beyond the plastic region. This lgh 
problem of trying to equate the impactor kinetic energy with the energy necessary to excavate 

the crater, as mentioned before. 

3.2 Cratering in TFE Teflon 

As was described above for cratering in aluminum, the CTH code has been u- d j° inv«hgate 
the cratering responses of a Teflon (TFE) target. Parameters varied have included impact 
velocity, projectile density, and target yield strength. In all cases the target was ^cjflytiuck 
to behave as an infinite body. For all cases the projectile was an aluminum sppere of diameter 
100 microns. As with the aluminum runs, calculations were done at 1 5 10 a»d 15 Vm/s, wtii 
yield values of 0.1Y„ Y„ and 10Y, (where Y t is the normal value for the Teflon target) and for 
densities of Pp /3, P p , and 3 Pp (where p p is the normal value of the alum num projectile). 
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The properties used for Teflon were: 


density * 2.17 g/cm 3 , bulk sound speed = 1.29 km/s, yield strength * 300 bars, s = 0.795, 
Poisson ratio = 0.433, and spall strength = 200 bars. The melt temperature is 327°C and melt 
energy is 82 cal/g. These are POD's "best guesses" for the material parameters, since the 
literature (Harper, 1992; Moses, 1978; Dean, 1992; Rice, 1980) reveals some inconsistences 
(some references even fail to distinguish between TFE and FEP Teflon, which have somewhat 
different properties). 

As with aluminum, the resulting data have been mapped and compared to equations of the form 

P/d, - constant (p/pJ A (p/YJ® u, c ( 96 ) 

and djd, - constant tp/p J/(l +(p (p/YJ" «, c (c/c^ D (97) 

Figures 33 to 42 show the data in log-log form. It is immediately clear that the indexes are not 
constants and that they vary with the impact speed, as occurred for aluminum. This effect is 
strongest at the lowest speed of 1 km/s (which normally would not be considered as 
hypervelocity) and slowly asymptotes to steadier values at the higher impact speeds. 

To date, the best overall fits for Teflon TFE are given by: 

for P/d,: constant - 0.68, A *= 0.49, B - 0.23, C - 0.475 

(Note that these indexes are similar to those in Equation 28, except for the density index and the 
leading constant). 

for dj/d,: constant - 0.18, A - 0.333, B = 0.20, C - 0.58, D - 0.333 

These fits are primarily for data above 5 km/s impact speeds. 

Of particular interest are the plots of (P/d,) and (d^d,) versus projectile density and speed 
(Figures 37 and 38). Note how for die crater diameters the plots "curve over" as projectile 
density increases. This is the form predicted by Equation 21. For the penetration depths the 
plots follow a more direct power law, which is as predicted by Equation 28. Figures 41 and 42 
show the ratios P/d e . 

Figures 43 to 51 show some of the geometries of the crater for the 10 km/s impact velocity case. 
The dark spots on the impact axis are Lagrangian points, as described below. Note the dynamic 
variations in die shape of the crater and the repeated sequence of lip spallation. The final crater 
shape has a surface spall region surrounding the crater proper. This is very similar to the features 
observed experimentally, as shown in Figure 4. 
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A problem was met when ttying to detemun. fcect^r region 

For some of the larger cmters the J—M-- - Conseqtmndy. more than 

displaying re-entrMit shapes due to bp apritoon ™ in the initial target surface 

C ^ f A. PTH runs reveals many interesting facts. It is observed that, as for the aluminum 

Scrutiny of the CTH occurs (in Fi 49 ) before the maximum crater diameter 

f 85 ®’ e S^te^ure 51) and that the ratio depth/diameter is not 0.5 in general {?.«., die crater 
has formed (in Figure 5 1 Umd ~ ^ ^ indica ted a "central ’pip*" at the 

^'^fthe°^.”^ith tfte'deepest pint of the crater slighUy off the impact «ds). Various 

ST ^(SSSi PO-)- Placed " 

r^a p ;« 

^ted^ve 2 ^ Llgrmgial points ^re LI through L7 at initial depths of 0.0 (unpact surface 

a m ^™ d **°* 500 microns gi™g p /d^ 
diameter was 720 microns. Thus d/d, - 7.2, and P/d* - «•«"• m an c^so 
depth and diameter are relative to the initial target surface. 

Immediately upon impact the peak suets is 

^Sss^SSSss^asssss: 

asymptotes towards a 1/r 24 law with increasing distance, 
fundamental concept in POD's derivation of its scaling laws. 

The induced compressive suai» during dte 1-D comhtionm Jj 

JWL increased £ di redo W/(W-u) 
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which is 6.8. The plots indicate a veiy small elastic recoveiy. A stress gradient is attached to 
the rear of the shock pulse at a value of about the target yield strength. 


As with the aluminum runs the CTH calculations demonstrate that die yield strength of die target 
is indeed a major player in the crater formation. Further, a large region (extending to beyond 
twice the crater depth) of the target is put into permanent plastic strain. This effect is most 
clearly seen for the deepest Lagrangian points (Figures 58 and 59) where it is seen that the 
material is ultimately left with a permanent compression of about the yield stress. This is the 
response expected for a material which has been subjected to a pseudo 1-D stress wave which 
mtceeds the Hugoniot Elastic Limit (HEL). Note that HEL = Y(l-u)/(l-2u) where u is the 
Poisam ratio. For Teflon u - 0.433, thus HEL = 4.23Y (i.e., 1.27 kbars for Teflon TFE). Upon 
shocking, the material returns to a state of zero stress in the direction of the stress wave but has 
lateral compressive stresses equal to die yield stress. Thus the pressure becomes 2Y/3 At these 
large radii the hoop stress is very small and provides only a small correction. 


A considerable amount of strain energy is stored within this plastic region beyond the crater, and 
energy also exists in the propagating elastic wave beyond the plastic region. As stated above 

this highlights the problem of trying to equate the impactor kinetic energy with the enerev’ 
necessary to excavate the crater. ** 


3*3 Perforations in Alumin um 


CTH has been used to investigate the perforation responses of an aluminum 6061-T6 target 
Panuneters varied have included impact velocity and the target thickness. For all cases the 
projejmle was a sphere of diameter 100 microns. As with the aluminum cratering runs 

calculations were done at 1, 5, 10 and 15 km/s, with the normal yield value of Y„ and normal 
densities of p p and 

The properties used for A1 6061 -T6 were: 


d *? sily „ g/C ^’ bulk 501111(1 speed = 51 km/s - y»*M strength = 2.7 kbars, s = 1.4, Poisson 

ratio * 0.33, spall strength = 3.1 kbars. 

The resulting data have been mapped and compared to equations of the form 

T/d p - const (p p /p t ) A (p/Y,) B u, c (98) 

and to POD's Equation 59. 


Figures 60 and 61 show the data in both linear and log-log form. It is immediately clear that the 
velocity index is not constant and slowly varies with the impact speed. This variation is strongest 
at the highest speeds above 10 km/s. 


To date, the best overall fits to Equation 98 for aluminum are given by: 
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constant « 0.133, A - 1/3, B - 1/3, C - 0.75 


The direct application of Equation 59 is shown in Figures 60 and 61. The responses are similar 
to the expectations. However, locating die Ballistic Limit threshold requires aiany calculation? 
(e.g., at least two for a given impact speed, one above and one below the thicshold). For this 
reason POD has yet to fully investigate the influence of density, yield strength *«id spall strength 
However, we anticipate responses similar to those observed for cratering. 

Figures 62 to 67 show some of die geometries of the crater for die 10 kni/s imp act velocity case, 
for a case of "not quite" causing perforation. The dark spots on the impact axis are Lagrangian 
points, as described below. These figures are split, with the right side showing geometry only 
and die left side giving stress contours as per the bar legend to the right of each figure. Note the 
dynamic variations in die shape of the crater and die repeated sequence of rear s urface spallation. 

The CTH runs reveal many interesting facts. For the non-perforation case (300 micron thick 
target) it is observed that the maximum depth of the crater occurs (in Figure 65) before the 
maximum crater diameter has formed, and that the ratio depth/diameter is not 0.5 in general (i.e., 
die crater is generally not hemispherical). Various tracer points (Lagrangian pc alts) were placed 
along the impact axis. For these points the time histories of pressure, velocity aad position were 
obtained. Figures 68 to 74 show such data for the case of a 100 micron A1 ball hitting a 500 
micron thick A1 target at 10 km/s. The Lagrangian points were LI through L7 at initial depths 
of 0.0 (impact surface of die projectile), 0.0 (surface of die target), 110 microns, 220 microns, 
340 microns, 460 microns and 570 microns, respectively. Of these, the points LI to L5 were 
within the finally formed crater which developed a maximum depth of 350 micions, giving P/dp 
= 3.5. Note that this depth is larger than the corresponding thick target case, where the depth 
was 300 microns. Also note the thin wall remaining at the bottom of the crater. The final crater 
diameter was 500 microns, which is identical to die thick target case. Thus &Jc = 5.0, and P/d e 
= 0.70. In all cases the definitions of depth and diameter are relative to the initial target surface. 

Immediately upon impact the peak stress is about 900 kbars. This is lower than use expected 1-D 
stress of 1.63 Mbars (perhaps owing to the artificial viscosity in die code?) but is considerably 
higher than the corresponding Bernoulli stress of 337.5 kbars. This peak stress rapidly drops 
through the Bernoulli value to a very low stress of the order of die plastic yield value, whereupon 
the stress drops much more slowly. Some residual stress reverberation ringing can be seen on 
die plots. The initial shock wave speed is about 12 km/s in the al uminum at die 1-D stress, while 
die release wave speeds are even faster. Thus the 1-D stress lasts for less than 10 ns and 
immediately collapses toward the Bernoulli stress. The peak stress then propagates at a wave 
speed of 6 km/s which is asymptoting to the dilatational speed for aluminum, and decays in value 
with distance. This decay roughly obeys a 1/r law near the impact point but steadily asymptotes 
towards a 1/r 2 law with increasing distance. 

The induced compressive strain during the 1-D condition is given by (u/W), where u is the 
induced particle speed change and W is the shock wave speed. Thus for a speed u of about 5 
km/s and a W of 12 km/s, the strain is 0.417 and the density has increased by the ratio W/(W-u), 
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which is 1.7. 
for this state. 


Thus the concept of treating the materials as "incompressible" is clearly not correct 
Thus the conventional Bernoulli logic applies only after the shock moves away 


from the impact site. 


The perforation case (400 micron thick target) is shown in Figures 75 to 80 (geometiy versus 
time), and in Figures 81 to 89 (Lagrangian tracer points). Hie figures clearly show that die crater 
mouth forms identically to that for the case of no perforation. The major difference is the deeper 
spall leading to crater break through. Note how the bottom of the crater has a conical shape. 
These effects have been observed experimentally for impacts which just cause perforations. 


3.4 Perforations in TFE Teflon 

CTH has also been used to investigate the perforation responses of a TFE Teflon target. 
Parameters varied have included impact velocity and the target thickness. For all cases the 
projectile was a sphere of diameter 100 microns. As with the Teflon cratering runs calculations 
were done at 1, 5, 10 and 15 km/s, with the normal yield value of Y t and normal density of p„ 
and p,. The properties used for A1 6061-T6 and TFE Teflon were the same as for the other 

calculations. 


The resulting data have been mapped and compared to equations of the form 
T/dp * const (p/p,) A (p/YJ B u, c 


(99) 


and to POD's Equation 59. 

The direct applications of Equation 59, together with the CTH data in both linear and log-log 
form are shown in Figures 90, 91 and 92. It is immediately dear that the velocity index is not 
constant and slowly varies with the impact speed. This variation is strongest at the highest 
speeds above 10 km/s. 

To the best overall fits to Equation 99 for TFE Teflon are given by. 


constant m 0.25, A = 1/3, B ** 1/3, C * 0.71 

Figure 90 shows a comparison of the CTH results against POD's Equation 59 for N = 2. As 
expected and explained above, the equation overpredicts the CTH data for this case. The 
explanation above stated that N = 2.4 should be the correct case for TFE Teflon and Figures 
and 92 confirm this. Locating the Ballistic Limit threshold requires many calculattons e.g at 
least two for a given impact speed, one above and one below the threshold), consequently PO 
has yet to fully investigate the influence of density, yield strength and spall strength. However, 
we anticipate responses similar to those observed for cratering. 


For the non-perforation case it is again observed that the maximum depth of the crater occurs 
before the maximum crater diameter has formed, that the ratio depth/diameter is not 0.5 in 
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general, and that in general die crater depth is larger than die corresponding t.ack target case. 
In addition, we observed dynamic variations in the shape of the crater and repeated sequences 
of lip spallation, and die CTH runs sometimes indicated a "central 'pip" at die bottom of the 
craters with the deepest part of die crater slightiy off the impact axis. The final crater shape has 
a surface spall region surrounding the crater proper. This is very similar to the features observed 
experimentally, as shown in Figure 4. In all cases die definitions of depth t ad diameter are 
relative to the initial target surface. 

Figures 93 to 101 show some of the geometries for the 10 km/s impact velouty case for the 
perforation case (700 micron thick target). The dark spots on die impact axi: are Lagrangian 
points. Note the dynamic variations in the shape of die crater and the repeated sequence of rear 
surface spallation. Also note the odd crater shape with local radial cracks, very similar to those 
observed by HOrz (1992) and shown in Figure 4. 

The crater mouth forms identically to that for die case of no perforation. The ajor difference 
is the deeper spall leading to crater break through. Note how the bottom of die crater has a 
conical shape. These effects have been observed experimentally for impacts which just cause 
perforations. 

Various tracer points (Lagrangian points) were placed along the impact axis. For these points 
the time histories of pressure, velocity and position were obtained. The I agran; : an points were 
LI through L9 at initial depths of 0.0 (impact surface of the projectile), 0.0 ^surface of the 
target), 110 microns, 220 microns, 340 microns, 460 microns, 570 microns, 680 microns and 800 
microns, respectively. Figures 102 to 108 show die data for Lagrangian points LI through L7 
for die case of a 100 micron A1 ball hitting a 700 micron thick TFE Teflon target at 10 km/s. 
Note from these Lagrangian points how slowly the perforation occurs with a very long-lived 
gradual downstream motion of the crater bottom. 

3.5 Oblique Impacts into Aluminum 

The CTH code has been used to investigate die oblique impact responses of an < uminum 6061- 
T6 target. Because such calculations are 3-D (and hence require much computer memory and 
run time) the only parameter varied was the angle of incidence, which was set at 70° and 50°. 
The projectile was a sphere of diameter 0.5 cm, and die calculation was done at .5 km/s, with the 
normal properties used for A1 6061 -T6. 

Figures 109 to 115 show the 70° impact as an isometric view. The projectile approaches the 
target from die top left, proceeds to cause a crater, and then ricochets off to the right. Figures 
116 to 122 show the same impact but in cross-section in die plane of die impat or (r e., a view 
"sideways" at the collision), while Figure 123 shows the orthogonal cro:s-sect: n (i.e., a view 
from "behind" the projectile). 

The figures clearly show how the crater develops a significant asymmetry, especially on the 
downstream side. The projectile is seen to "ride up" the extended crater wal and to slowly 
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stretch and fragment. For this impact, POD’s Equations 87 and 88 predict that the upper portion 
of the projectile will suffer a partial ricochet. The figures clearly show that, as predicted, die top 
of the projectile never actually enters the crater, but instead passes over the crater lips to impact 
die target further downstream. Note also how this portion of the projectile essentially moves 
parallel to the target surface. 

From die CTH data, die maximum crater depth is 0.42 cm. Thus P/dp = 0.8. The crater diameter 
measured in the orthogonal cross-section is d, * 1.0 cm. Thus d^dp = 2.0. For an impact speed 
of 5 km/s at 70° angle of incidence die component normal to die surface is 1.71 km/s. For this 
speed POD's Equations 21 and 28 give d t /d p « 2.0 and P/dp - 1.0, respectively. (Note that 
although POD's prediction for P is about 20% too high, no direct CTH calculations were done 
for an impact speed of 1.71 km/s normal to die surface. Thus die "exact" value expected is 
slighdy uncertain). 

For die 50° impact die prediction given before (Figure 10) suggests that die projectile should be 
on die limit of ricocheting. The CTH results confirm this. Figures 124 through 128 show a 
"side" view of the impact, and Figure 129 shows die "rear" view. It is seen that die top of the 
projectile initially just escapes its own crater at early time. However, the crater is still growing 
and succeeds in just catching the projectile fragments at later time, whereafter die projectile never 
subsequently escapes its own crater, although there is a lot of lip formation and ejecta The final 
lateral crater diameter is 1.44 cm (d,/dp = 2.88) while die depth is 0.634 cm (P/dp c 1.27). These 
values should be compared to 2.9 and l.S, respectively, which are the expectations for this case 
(normal impact speed component of 3.21 km/s). Thus die CTH calculations confirm that oblique 
cratering does obey die "cosine law" for depth of penetration and orthogonal crater diameter. 

An important point to note is that these particular CTH calculations were done for a much larger 
projectile (0.5 cm) than were all the previous calculations (100 microns) for normal impacts. Yet 
the scaling rules still apply! This demonstrates that (as POD previously stated) the 
"supralinearity" effect is not associated with hydrodynamics per se. However, this fact also 
suggests that the effect is not associated with changes in strain rates either. POD therefore 
believes that die Petch Law is still the most probable explanation for supralinearity. 

An issue raised in Section 1.4 concerned die possibility diat die cosine law breaks down for very 
high speed impacts, since die latter might be dominated by energy rather than momentum. 
However, none of POD’s analyses invoke energy as a criterion for cratering. Based on this fact, 
POD presendy believes that the cosine law continues to apply even for the highest impact speeds. 
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4.0 COMPARISONS WITH OTHER SCALING LAWS 

4.1 Cratering in Aluminum 

A large number of existing scaling laws describing either penetration depth and/c r crater diameter 
have been identified and compared. Christiansen (1992a,b) has given quotes for the Cour-Palais 
equations. Schonberg (1989) has also listed many of these laws. For convenience we include 
die equations from his report, together with his list of references, in Appencux B (Note dial 
Schonberg's data accidendy misquotes the Cour-Palais equations, giving a velocity index of 4/3 
instead of 2/3. For this reason POD uses die Christiansen quotes only). Figui# 130 compares 
die results of POD's Equations 21 and 28 with CTH calculations. We see that for a projectile 
and target both of A1 6061 -T6 the equations give good fits to the CTH predictions. Figures 13 1 
to 141 show comparisons between POD's equations and these other scaling laws for the same 
impact conditions. 

Of all the existing equations the one by Sedgwick et al. (1978) seems to be tho closest to the 
CTH data for both penetration depth and crater diameter. The good agreement with Sedgwick's 
equations is not entirely unexpected, since his equation is itself based on fits to calculations done 
with die HELP hydrodynamic code. If anything, this merely proves that CTH nd HELP give 
similar answers (both correct or both incorrect?). 

We see that the scaling laws of Cour-Palais (1985) (for dp = 1 cm) and Bruce 1962) are very 
close to POD's equations. Cour-Palais' equation utilizes a supralinearity term c 30S6 , and is 

P/dp = 5.24 dp 00 * (Pp/p / 5 (l/Hf 23 (uocose/coj 2 * (100) 

where H is die Brinell hardness number. 

It is of interest to note the good agreement of his equation for unit diameter. As Mated in Section 
2.3.9, POD believes the supralinear index effectively goes to zero for a 1 cm diaiaeter projectile. 

Unfortunately, the later data of Bruce (1979) gives a much poorer fit Other reasonable fits are 
those by Dunn, Goodman-Liles, Sawle, and (for crater diameter only) Summer#. The scaling 
laws of Christman, Herrmann-Jones, Sorenson, and Summers-Charters give much poorer fits. 
Disconcertingly, the law of Summers-Charters (P/dp) shows die wrong power index versus impact 
speed, indicating an index greater than unity. This is contraiy to all other law=, experimental 
observations and physical logic, except at very low impact speeds. It is possible that this is 
another misquote by Schonberg. POD has not yet checked this possibility. 

With regard to experimental data there are some ambiguities since not all resear-ners have used 
consistent definitions of crater diameter and depth. This would help explain why differing "fits" 
to similar experiments are often seen to be "off-set" from one another on the plots. The most 
respected definitions (Cour-Palais) are those that reference die original target surface, and POD's 
analysis and measurements from the CTH calculations have used these definitions. 
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4.2 Cratering in TFE Teflon 

The results of POD's Equations 21 and 28 were shown in Figure 52. We see that for a projectile 
of A1 6061-T6 these equations give good fits to the CTH predictions, being about 18% low 
across die entire velocity range. A large number of existing scaling laws describing either 
penetration depth and/or crater diameter have been identified and compared. Schonberg (1989) 
has listed many of these laws. For convenience, we list die appropriate laws, and Schonberg's 
references in Appendix B. However, these laws are mostly based on observation for aluminum 
targets. Therefore, those equations which do not explicidy contain material yield strengths cannot 
be used to describe Teflon targets. Consequendy, only those laws which do contain the yield 
paramf*t* r have been used. Figures 142 through 146 show die predictions from these other 
scaling laws for die case of an AL 6061-T6 projectile into Teflon. 

If these other scaling laws were truly "generic" they should fit the Teflon data just as well as they 
did die aluminum data, when die appropriate material properties are used. It is clear, however, 
that this does not happen. POD's equations, however, do gives good fits to both Teflon and 
aluminum data. 

Only the equations of Sedgwick, Bruce (his 1979 version), Dunn and Sorenson can be used for 
these comparisons, since none of the others contain the yield parameter (although some contain 
hardness numbers). 

The Cour-Palais equation makes use of Brinell Hardness of die target, rather than the yield 
strength. This is unfortunate, since the Cour-Palais equation is frequently used by NASA for 
cratering predictions. The problem with hardness numbers is that they do not direcdy relate to 
other material properties, such as yield, and are difficult to compare if different scales are used 
for different materials. Thus while aluminum has a well quoted Brinell hardness number (BHN 
95 for 6061 -T6), TFE Teflon has a Shore hardness of D52. An approximate equivalent is a 
Brinell value of about BHN5 (or Rockwell R60) (Harper, 1992). Using the Cour-Palais equation 
with this hardness value, with a Teflon density of 2.17 g/cm 3 , and a sound speed of 1.29 km/s, 
die predictions are too large for penetration depths or diameters, as seen in Figure 143. 

The equation by Sedgwick et al. (1978) is reasonably close to die CTH data for penetration 
d e pth , but it significandy overpredicts die crater diameter. Sedgwick's equations give good fits 
for al uminum, but not for Teflon. Since his equation is itself based on fits to calculations done 
with die HELP hydrodynamic code for aluminum, this demonstrates that die latter fit is somewhat 
fortuitous. POD believes this is an example of wrongly choosing the "pi" groups for the 
equation, as discussed earlier. 

We see that the scaling laws of Bruce (1979) give veiy poor fits (grossly overpredicting), for 
either depth or diameter. The equations of Bruce (1979) give very poor fits to aluminum data, 
also. Dunn's equations also give gross overpredictions for both depth and diameter for Teflon, 
but give reasonable fits for aluminum. Lasdy, Sorenson’s (1962) equations give a good fit for 
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Teflon depth daht, but overpredtc, the diameter data Sorenson's equates give gross 
underpredictions for both depths and dimeters for alummum. 

■ggssaaagag- ,, "*‘ 

POD concludes that, unlike the POD equMmns presented here, practically none f the previous 
existing scaling laws can be considered generic. 

43 Perforations in Aluminum 

convenience, we list the app P > n f an A1 6061 -T6 projectile into A1 
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how well either of these equations fit the CTH data. However, both equations have the highest 
of all the various velocity indexes, suggesting over-predictions at higher speeds. 

For predictions of die perforation hole sizes, equations have been given by Maiden et al (1964) 
McHugh (1962), Sawle (1969), and Brown (1970). All of these equations have the property of 
predicting significandy ever larger holes as die impact velocity increases. Further, only that by 
Sawle includes material densities, most only use linear dimensions! 

Herrmann and Wilbeck (1986) point out that such equations rarely fit a wide range of data. The 
equations by Sawle, Maiden and Brown imply that die hole size increases with almost a (T/d.) M 
rule, but the data is not well confirmed. Brown observed "odd non-linearities" as the target 
thickness decreased to zero, but no rational explanation was given. 

POD s Equations 75 and 76 indicated a roughly simple linear rule for hole size versus (T/d.) 
when die latter was veiy small, but a reducing index as T/d„ increased. ^ 

4.4 Perforations in Teflon 

As described above for aluminum, several existing scaling laws describing perforation, 
specifically the Ballistic Limit, have been identified and compared. McDonnell and Sullivan 
(1992) have listed many of these laws. For convenience, we list the appropriate laws and 
McDonnell's references in Appendix B. Figures 149 and 150 show the predictions from the’laws 
for the case of an Al 6061-T6 projectile into TFE Teflon. Also shown are the results of POD's 
Equation 59. We have previously shown in Figures 91 and 92 that, for a projectile of Al 6061- 
T6, this equation gives a good fit to the CTH predictions. 

In Figure 149, it is seen that of the various equations by McDonnell (1992), his Equations 6, 7, 
and 10 all fit die CTH data (slightiy low) for Al 6061-T6 projectile and TFE Teflon target! 
provided the value of 800 bars is used for o A1 (as done by McDonnell), q is set to 0.2 kbars (as 
with the CTH data), and dp = 1.0 cm (to ••remove" the supralinear term, as suggested by POD 
m Section 2.3.9 for such large projectiles). However, his Equation 11 McDonnell, 1979) is a 
poor fit, tending to give very low predictions. McDonnell assumes that the tensile strength 
should be used, whereas POD concludes that both yield and tensile strengths are needed, as per 
Equation 59. All three of these McDonnell’s equations are close fits to die CTH data. However, 
since McDonnell recommends (in his report) the use of his Equation 10, and since it fit POD's 
Equation 59 the best, we have used this equation in Figure 150 to compare with the predictions 
of scaling laws from other investigators. 

The equation by Naumann (1966) overpredicts die CTH data by a factor of about two, and is 
clearly not good for aluminum into TFE Teflon. In Section 4.3 it was shown that Naumann was 
also high by a factor of about three for aluminum into aluminum. 

The equation by Cour-Palais (1979) is too low by a factor of about 0.6 for an Al/Teflon impact. 
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Both the equations of Fish and Summers (1965) and of Pailer and GrOn (1980) involve a term, 
e , describing ductility. Since this quantity is not well defined it is difficult to determine how 
well either of these equations fit the CTH data. For TFE Teflon, we took this quantity as an 
average of the values quoted by Shackelford and Alexander (1992). Thus we used 6= 275 %. 
However, for the Al/Teflon impact. Fish and Summers predict low (similar to Ccyr-Palais), while 
Pailer and GrQn are very close to POD's Equation 59. 

If these other scaling laws were truly "generic” they should fit the Teflon data jut, as well as they 
did die aluminu m data, when the appropriate material properties are used. It is clear, however, 
that this does not happen. POD's equation, however, does give good fits to both Teflon and 
aluminum data. 
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5.0 COMPARISONS WITH EXPERIMENTAL DATA 

5.1 Cratering in Aluminum 

Some recent well character! zed experimental impact data have been provided by HOrz (1992). 
These experiments have made use of a gas gun at NASA Johnson Space Center to accelerate 
projectiles of soda-lime glass at aluminum targets of A1 1100 and A1 6061-T6. The impact 
speeds were all about 5. 8-5.9 km/s and were normal to the target surfaces. The projectiles were 
of 50, 150, 1000 and 3175 microns diameter. H6rz carefully measured the crater diameters for 
each experiment, using the initial target surface as die reference plane. The primary purpose of 
these experiments was to study perforations, consequently only a small fraction of the data were 
for craters in thick targets. 

For the thick A1 1100 targets die values of (d/dp) versus projectile diameter were found to be 
4.33 (3175 microns), 5.01 (1000 microns), 4.92 (150 microns) and 3.2 (50 microns). For the A1 
6061-T6 the ratio was 3.35 (3175 microns). This latter value is about 15% lower than the POD 
CTH calculation for an AI/A1 impact. Soda-lime glass has a density of about 2.2 g/cm 3 and a 
slightly lower sound speed than aluminum. Accordingly, Equation 21 predicts only a very small 
difference (i.e., a few percent) versus the Al/Al impact. 

Although the d ata base is small and therefore subject to some error, a distinct drop in crater 
diameter is seen for the smallest projectiles versus the larger ones. The larger projectiles tend 
to give a systematic trend in crater diameter, except for the largest ones. Comparing the 50 
micron projectile data with the 3175 micron data there is an apparent supralinear index of about 
0.073, while comparing the 50 micron data with die 1000 micron data gives an index of 0.149. 
These apparent indexes are both larger than die Cour-Palais quote of 0.056 and may be a 
consequence of die small data base. 

Comparing the results for A1 1 100 versus those for A1 6061-T6 we see that the crater diameters 
were in the ratio of 4.33 to 3.35, i.e., 1.293. The only significant difference in these two 
al uminums are the yield strengths. Equation 21 gives an index of -0.2857 for crater diameter 
versus yield strength, while die fit to die CTH data gave an index of -0.258. Thus to explain die 
experimental data, assuming the yield strength of A1 6061-T6 is 2.7 kbars, requires the A1 1100 
to have a yield strength of either 1.098 or 0.997 kbars. Engineering data (Shackelford, 1992) 
identifies this alu minum to be of temper H12 (i.e., A1 1100 H12), which is one of the softer 
aluminums. 

Recent cratering data have also been provided by Christiansen (1992c) for aluminum projectiles 
into A1 6061-T6. This work was done at the NASA JSC HTT-F gas gun facility, as part of a 
careful study of oblique impact responses. For a normal impact with an impact speed of 6.83 
km/s the experiments give P/dp = 2.1 and d/dp = 4.55. These values are close to both the CTH 
predictions and POD's predictions using Equations 28 and 21, respectively. CTH gives P/dp = 

2.2 and d/dp = 4.0, while POD's equations give P/dp = 2.2 and d/dp = 4.1. 
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5.2 Cratering in Teflon 

There is only a limited data base for impacts into Teflon TFE. Some recent weil characterized 
experimental impact data have been provided by H6rz (1992) for project les of oda-lime glass 
and Teflon TFE targets. The impact speeds were all about 5. 8-6. 3 km/s and were normal to the 
target surfaces. The projectiles were of 50, 150, 1000 and 3175 microns diameter, but only the 
largest ones were used for cratering in thick targets. Hflrz carefully measured the crater diameters 
for each experiment, using the initial target surface as the reference plane. The primary purpose 
of these experiments was to study perforations, consequently only a small fraction of die data 
were for craters in thick targets. H6rz also sent photographs to POD of cross-sections through 
the Teflon samples, shown in Figure 4. It is evident from these pictures that the c aters in Teflon 
are not smooth-surfaced as with a ductile metal target. Instead, the craters suffer from rough 
surfaces and show signs of surface fracture, some radial cracking, and also have a "hairy" 
morphology indicating "strings" of partially melted/resolidified material ( with the exception of 
the "strings" these morphological features are similar to the CTH example given above). Because 
of these features the "exact" value of crater diameter is difficult to define. For example, for the 
case of D/T f = 0 125 in Figure 4, H6rz quotes a crater diameter of 1.18 cm, for a 3175 micron 
soda-lime projectile into a one inch thick Teflon target (*.«., (d/d,,) = 3.72). POa's independent 
estimate of this diameter is 1.3 to 1.5 cm (i.e., (d/dp) = 4.09 to 4.72). based solely on the 
photograph. Thus a possible 27% error in die quote exists. This is consistent witi Harz's (1992) 
findings that the "standard" crater in Teflon was difficult to measure. 

Thus for die thick Teflon targets the values of (d/dp) were found to be in the range 3.72 to 4.72 
(3175 micron projectile). The CTH value is about 5.7 for an impact speed of $.3 km/s and a 
Teflon yield strength of 300 bars. This CTH value, for an Al/Teflon impact, is higher than the 
experimental value(s), for soda-lime glass/Teflon impacts, by a factor of between 1.21 to 1.53. 
Soda-lime glass has a density of about 2.2 g/cm 3 and a slighdy lower sound speed than 
aluminum, so Equation 21 predicts a very small difference (i.e., a few percent) versus the 
Al/Teflon’ impact. To make the CTH calculations agree with the experiments it would be 
necessary to either increase the assumed yield strength, or decrease the assumed bulk sound speed 
of die Teflon. 

With die exception of Sedgwick's 1978 law (which is close to the CTH penetration depth 
predictions for Teflon TFE) and Sorenson's 1962 law (which is close to the Equation 28 
penetration depth predictions for Teflon TFE), the other scaling laws (presented in Section 4.0) 
substantially overpredict the CTH values and, thus, the experimental data However, as 
previously stated, the Equation 21 and 28 predictions of crater diameter and penetration depth 
are consistendy 18% below the CTH values. Thus these equations predict values between die 
CTH values and Hdrz's experimental data. As shown in Figure 52, for the experimental velocity 
regime over which Harz found the Teflon (d/dp) values to be in the range of 3.72 to 4.72 
Equation 21 predicts (d/dp) values in the range of 4.6 to 4.95. Picking the highest vidue (at 6.3 
km/s) Equation 21 is higher than the experimental value(s) by a factor of between 1 .05 and 1.33, 
for an Al/Teflon impact. Since the soda-lime glass used in Harz's experiments has a lower 
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density and sound speed than aluminum. Equation 21 predicts a very small difference (i.e., a few 

percent) versus the Al/Teflon impact. 

Since only one projectile size was used there is no data concerning supralinearity. 

S3 Perforations in Aluminum 

Some recent experimental impact data have been provided by H6rz (1992) for projectiles of soda- 
lime glass and thin aluminum targets. The impact speeds were all about 5. 8-6.3 lcm/s and were 
normal to die target surfaces. The projectiles were of 50, 150, 1000 and 3175 microns diameter. 
However, only the largest projectiles were used to obtain data on A1 6061-T6 (modelled in this 
report). All projectile sizes were used for A1 1100 targets, however. H6rz carefully measured 
the crater diameters for each experiment, using die initial target surface as the reference plane. 
The primary purpose of these experiments was to study perforations. Harz also sent photographs 
to POD of front and rear views of the aluminum samples, shown in Figures 151 and 152. 

Although Harz used soda-lime glass instead of aluminum for the projectiles, POD's eq ua tions 
indicate that the difference in perforation responses is only a few percent. For A1 6061-T6 the 
experiments indicate that the Ballistic Limit occurs for T/dp of about 3.0. As shown in Figures 
60 and 61, this number is very close to the CTH predictions of 2.8 (perforation) to 3.3 (no 
perforation) for an impact speed of 6 km/s. Likewise, POD's Equation 59 gives a similar, but 
slighdy lower, value of 2.6. 

For A1 1 100 the major material properties which differ from A1 6061 -T6 are the yield and tensile 
strengths. Based on the cratering data we deduced that the temper of die metal was HI 2. 
Standard handbook quotes (Shackelford and Alexander, 1992) give die yield strength of A1 1 100 
H12 as 1.0 kbar and the tensile strength as 1.1 kbars, versus 2.7 kbars and 3.1 kbars, respectively, 
for 6061-T6. Equation 59 suggests that this will increase the limiting value of T/dp to about 3.5 
at 6 km/s. This value is consistent with the experimental data 

Comparing Hdrz's data as a function of projectile diameter reveals that a supralinear behavior 
exists. Comparing die 50 micron results versus die 3175 micron results, the data indicates an 
index of about 0.07, which is a little higher than the Cour-Palais quote of 0.056, and both values 
are consistent with die supralinearity indexes expected from POD's approach in Section 2.3.9. 

5.4 Perforations in Teflon 

Recent experimental impact data have also been provided by HOrz (1992) for projectiles of soda- 
lime glass and thin TFE Teflon targets. The impact speeds were all about 5.8-6.3 km/s and were 
normal to die target surfaces. The projectiles were of 150, 1000 and 3175 microns diameter. 
H6rz carefully measured the crater diameters for each experiment, using the initial target surface 
as the reference plane. The primary purpose of these experiments was to study perforations. 
Harz also sent a photograph to POD of cross-sectional views of die Teflon samples, shown in 
Figure 4. 


56 



Although HOrz used soda-lime glass instead of aluminum for the projectiles, J>U>s 
indicate that the difference in perforation responses is only a few percent. For TTE Teflon the 
experiments indicate that the Ballistic Limit occurs for T/d, of about SX* 0.25 AiAmib 
F igures 91 and 92, this number is veiy close to the CTH predictions of 5.3 (perforation) to 6.4 
(mfperforation) for an impact speed of 6 km/s. Likewise, POD’s Equation 59 eves a similar 
but slightly lower, value of 5.1. All of the other investigators’ (e.g., Mc^nneli Cour-P^s) 
predict much lower values (the closest being Pailer and GrOn with a value of about 4.7), with the 
exception of Naumann (who substantially overpredicts a value of about 7.9). 

Harz's data as a function of projectile diameter is presently insufficient to allow determination 
of supralinear behavior. 

5.5 Oblique Impacts in Thick Aluminum 

Some recent careful studies of oblique impacts of aluminum into duck A1 5061- o targets have 
been done by Christiansen (1992c) at the NASA JSC HIT-F gas gun facility Christiansen 
studied the responses over an obliquity range of from 0 incidence to 88 incidence, using sm 
incremental angular changes at the high obliquities. The experiments were all done w.A an 
impact speed of about 6.75 ± 0.2 km/s. The data show that the value of P/dp roughly obeys a 
(cose) 2 * law, as expected, but tends to drop somewhat faster at very large angles of incidence 
above 70°. This behavior is similar to POD's prediction of Equation 83 which is s own m Figure 
8 The data also show that the value of d./dp roughly obeys a (cos6) * law, Mis an even 
better fit to a (cosG) 057U law as given by POD's Equation 82 and shown m Figure 8. 

Of great interest is the data for <Wd^ . where is the attended thameierl ^ 

L is ,he perpendicular diameter. For angles of incidence below about 69 this ratio is 
essentially 1 0 indicating that the craters are basically axisymmetric. However, for larger angles 
the ratio rapidly increases. To a first approximation the ratio fie *e rule O.56nm0 over Ae 
angular ra nge 60° to 88°. This behavior should be compared with die ian0 -ogic of PO 

89 which describes die ability of the top of the projectile to just escape its own crater 
An interpretation is that by just failing to escape its own crater the projectile is causing the 
downstream stretching of the crater. 

5.6 Oblique Impacts in Thin Aluminum 

Unfortunately, there is little data concerning the definition of the Ballistic Limit for dun target 
perforations L functions of obliquity. Instead, most experimenterehave used 
targets to study the effects of obliquity on ricochet of the projectile and perforation d etas .eje cta. 
Christiansen (1992c) and Schonberg (1988, 1989) have done such studies. As with cratenng 
is observed that ricocheting occurs above about 60° angle of incidence. 
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6.0 OTHER MATERIAL COMPARISONS 

6.1 Copper Cratering 

POD has recently received details of work done by Wingate et al. (1992) at Los Alamos National 
Laboratory (LANL), presented at the 1992 HyperVelocity Impact Symposi um (HVIS). The work 
involves copper on copper impacts, and compares four code predictions. The codes are: EPIC, 
MESA, SPH and CALE, and experimental data is also compared. The following table lists die* 
codes' results and the POD predictions. 

The calculations are for an impact at 6 km/s. The projectile diameter was d_ * 0.4747 cm (0 5g) 
Properties for copper were: p p = p, = 8.93 g/cm 3 , c* = 3.94 km/s, s = 1.49, Y = 24 kbars 
To compute our values POD used Equations 28 for P and 21 for d,. 

RESULTS 


QUANTITY 

EXPERIMENT 

EPIC 

MESA 

SPH 

CALE 

POD 

P (CM) 

1.4 

1.8 

1.59 

1.73 

1.51 

1.55 

D c (CM) 

2.54 

2.4 

2.8 

2.6 

2.44 

2.71 

P/D c 

0.55 

0.75 

0.57 

0.67 

0.62 

0.572 


We observe that POD's predictions are close to the experimental data. Also note that the 
variations in the code answers are themselves about 19% (for P), 17% (for dj and 32% (for P/dA 
The ratios for the POD values versus experiment are: c 

1.107 (for P), 1.07 (for dj and 1.04 (for P/dJ. 

Part of Wingate's work was to explain supralinearity for small (micron size) projectiles. To do 
so he invoked strain-rate hardening and proposed that the effective yield strength of copper acted 
as if 5 times larger than normal, and thus was set at 12 kbars. This increased yield value reduced 
foe code predictions for crater volume by a factor of 4.1 (EPIC), 4.4 (MESA) and 3.3 (SPH) 
The POD prediction is 3.973 (Equation [21] 3 ) for the same higher yield. Note that LANL did 

not actually use a strain-rate model, they merely increased the yield value in die elastic-plastic 
model. 

6JL Lead Cratering 

As part of his recent studies of cratering and perforation, H6rz also used lead targets with soda- 
mie glass projectiles. Hdrz also sent photographs to POD of cross-sectional, front and rear views 
of the lead samples, shown in Figures 153, 154 and 155. Using Equation 21 for crater diameter 
we can find the ratio of lead crater diameter versus aluminum 6061-T6 crater diameter for 
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constant projectile size and impact speed. Assuming a lead density of 1 1 35 g/c -i 3 , sound speed 
of 2 krn/s, and yield strength of 0.13 kbars, we predict that lead craters are approximately 1.34 
times larger than aluminum craters. HOrz's data gives a ratio of about 1.48. T us Equation 21 
appears to give a good prediction. 
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7.0 CONCLUSIONS 


7.1 Cratering 

Based on die present work POD believes it has a strong insight into cratering laws. Clearly die 
concept of hemispherical craters is rarely correct, and it appears that diameter is governed by 
different rules from those giving depth, although the craters do asymptote towards hemispheres 
for die higher impact speeds. The work demonstrates that target strength is a strong driver for 
crater size. While increasing die ratio of projectile density to target density always increases 
crater sizes die responses are not simple power laws. Likewise die power index for impact speed 
is not the simple 2/3rd value. Overall, the data suggest that several of the indexes are "coupled" 
(e.g., die index for density is itself dependent on die velocity, etc.). 

POD believes it has identified the source of die supralinearity observed in cratering, and has 
shown the effect to be related to material strength behavior, in particular that of the target. The 
resulting analysis suggests that the supralinear effect is really a small-size downgrading and that 
die effect essentially vanishes for craters larger than about 1 cm diameter. Although insufficient 
material property data exists for Teflon TFE, POD believes the size effect to be roughly the same 
as for aluminum. Accordingly, projectiles of size 3175 micron (1/8 inch) are close to the limit 
where supralinearity asymptotes. Unfortunately, there were no quotes by HOrz for crater 
diameters for smaller projectiles. 

POD has considered the possible effects of momentum enhancement for grossly vaporizing 
projectiles (e.g., for very high impact speeds), and has concluded that this enhancement is 
generally only a few percent. Accordingly, the effect is minor and no obvious "step jumps" in 
cratering responses are expected as materials vaporize. The CTH aluminum and Teflon 
calculations do not indicate any "vapor" effects on either crater diameter or penetration depth up 
to 15 km/s impact speeds, even though partial vaporization of both die aluminum projectile and 
some of the aluminum and Teflon targets occurs at the higher speeds. (Note that both the Mie- 
Gruneisen and ANEOS equations of state were tried in the CTH runs with little difference in 
results). 

From the above efforts, POD believes that, for aluminum, die physics-based scaling laws 
presented in Equations 21 and 28 should be used for making predictions of, or interpreting data 
from, crater diameters and depths, respectively. An alternative would be to use Sedgwick's 1978 
scaling law for penetration depth. However, if Sedgwick's law is used, it must be borne in mind 
that Sedgwick does not account for the difference between crater diameter and penetration depth. 
A further alternative would be to use the two CTH fits which do differentiate between diameter 
and depth. However, both of these alternatives are based purely on a hydrodynamic approach. 
For these reasons, POD presently recommends that its Equations 21 and 28 be used. 

Based on the present work, these scaling laws apply to both aluminum and Teflon, and also seem 
to work well for copper and lead. Equations 21 and 28 are dimensionless; however, these 
equations were derived from physics and thus suffer none of the irregularities (caused by 
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improper selection of "pi" terms) which hamper the other existing scaling laws In comparisons 
to CTH hydrodynamic code calculations, Equations 21 and 28 closely matches the predictions 
for al uminum and was within 18% of the predictions for Teflon. In addition, these equations 
predict trends in crater diameters and penetration depths, with varying target and projectile 
material properties and impact velocities, which match the trends predicted by CTH. In 
comparison with experimental data (HOrz, 1992), Equations 21 and 28 closely matched the data 
for aluminum impacts and was within a factor of 1.0S to 1.33 for Teflon impacts. 

This work has also shown that the maximum depth of die craters usually o :urs before the 
maximum crater diameter has formed. 

While alternative scaling laws exist, none of die previous existing scaling laws c m be considered 
generic since they fit impact data either for aluminum or for Teflon, but not fb= both materials. 
In addition, these other scaling laws are only applicable over limited velocity and material 
property regimes (i.e., they do not correcdy predict trends in crater diameter, or penetration 
depth). For the above reasons, POD strongly recommends that, for both aluminum and Teflon 
TFE, die physics-based scaling laws presented in Equations 21 and 28 should be used for making 
predictions of, or interpreting data from, crater diameters and depths, respectively. 

72 Perforations 

Based on the present work, POD believes it has an insight into perforatioi laws, and has 
developed a genetically applicable scaling law (Equation 39) for predicting th» Ballistic Limit 
perforation conditions for both aluminum and Teflon. We conclude that the equation describing 
die Ballistic Limit has two parts, one relating to die crater depth, and one /elating to the 
reflection of the shock pulse from the target rear surface. The former term d pends on yield 
strength and has roughly a 2/3 index for impact speed, while die latter term demands on tensile 
strength and has roughly a unit index for impact speed. Thus the work demongrates that both 
target yield strength and tensile strength are strong drivers for determining per -rations. Since 
most other researchers use only a single equation term, it is not surprising that gsibiguities arise 
with regard to which material strength term to use, and with regard to die correct speed index. 

POD believes the supralinearity effect also applies to perforation. 

From the above efforts, POD believes that, for aluminum and Teflon, th physics-based 
perforation law presented in Equation 59 should be used for making predictions of the Ballistic 
Limit. An alternative would be to use any of the three referenced McDonnell and Sullivan 1992 
scaling laws (Equations 6, 7 and 10) for perforations in aluminum. There is no clearly obvious 
"best case" among these three McDonnell equations. McDonnell's Equation 10 gives the closest 
fit to POD's Equation 39, and all three of these McDonnell equations give close fits to the CTH 
data A further alternative would be to use the CTH fit. However, these McDujinell and CTH 
fit equations are not based on the physics approach (as taken by POD) and may not be directly 
applicable to other materials. 
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Equation 59 is dimensionless; however, this equation was derived from physics and thus suffers 
none of the irregularities (caused by improper selection of "pi" terms) which hamper the other 
«wri«rin g scaling laws. In comparisons to CTH hydrodynamic code calculations for perforations. 
Equation 59 closely matched die predictions for aluminum and was within 10% of the predictions 
for Teflon. In comparison with experimental data (HOrz, 1992), Equation 59 closely matched the 
Hntn for al um inum impacts and was within a factor of 1.09 to 1.19 for Teflon impacts. 

This work has shown that, for die Ballistic Limit case, the crater diameter is essentially the same 
as for die semi-infinite target case. However, this work has also developed physics-based 
equations (Equations 72, 76 and 77) which predict both this effect, and the effect observed 
experimentally in ultra-thin targets (foils) where the crater diameter asymptotes to die projectile 
diameter. 

While alternative scaling laws exist, none of the previous existing scaling laws can be considered 
generic since they fit impact data either for aluminum or for Teflon, but not for both materials. 
POD's equation, however, does give good fits for both Teflon and aluminum data. In addition, 
these other scaling laws are only applicable over limited velocity and material property regimes 
{i.e., they do not correcdy predict trends in crater diameter or penetration depth). For the above 
reasons, POD strongly recommends that, for both aluminum and TFE Teflon, die physics-based 
perforation scaling law presented in Equation 59 should be used for making predictions of the 
Ballistic Limit. 

7 3 Oblique Impacts 

POD believes it has confirmed that for oblique impact the component of impact velocity normal 
to die target correcdy describes the target responses. POD has implemented a cos0 correction 
in its Equations 21, 28 and 59 to account for this phenomenon, giving final corrected Equations 
82, 83 and 84. POD recommends that these equations be used for impacts at any angle of 
incidence into both aluminum and TFE Teflon targets. 

POD has also indicated the important factors which determine the process of ricocheting of the 
projectile and has developed and validated criteria (Equations 87 and 88) for predicting partial 
and complete projectile ricochets. 

F inall y it should be noted that all of POD's equations were derived from physics using the logic 
of die conservation of momentum versus conservation of eneigy. 

7.4 Implications for LDEF 

The primary implication for LDEF (and for other returned spacecraft materials) is that the Cour- 
Palais equations should be used with care for data interpretations, since while these equations 
give good fits for thick aluminum targets, they give much poorer fits for Teflon and for 
perforations. For the Ballistic Limit for thin foils, die McDonnell equations (McDonnell and 
Sullivan, 1992) are the best existing scaling laws and can be used (although we recommend using 
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POD's equations in the future). However, the McDonnell equations should not be used for cases 
where the target thickness is large compared to the particle diameter. It should be noted that this 
latter is largely die case for LDEFs TFE Teflon blankets. For these reasons, we recommend that 
the LDEF d ata be interpreted using the POD equations presented in this series =jf reports. In 
addition, the Solar Maximum Mission data should be reinterpreted for determining die meteoroid 
and debris particle impactor sizes for use in developing the environment models 
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8.0 RECOMMENDATIONS FOR FUTURE WORK 


For future work POD recommends that cratering studies be done for other ductile and britde 
target materials, including plastics (e g., Kevlar, Mylar, Kapton). Also there remains a need to 
study the effects of projectile shape, and the progression of responses as die shape changes from 
plate-like through spherical to rod-like. 

POD believes that most of die uncertainties in die appropriate indexes could be resolved by more 
detailed studies of die "inner responses" of the CTH code calculations (e.g., using many tracer 
points to track out the behavior of stress, velocity and motion). This should better determine the 
"N+l" index in Equation 21. Likewise, such detailed studies done for die lower impact speeds 
would allow a better understanding of the regime of small craters, where the assumption of 
hemispherical shock waves is far from valid. This is die region where the responses change from 
die Bernoulli flow state into the elastic-plastic, and can occur at speeds of greater than 3 km/s 
for the case of strong ceramics. 

A parameter that needs tracking is the strain-rate effect on cratering. POD's present analyses used 
only the standard elastic-plastic models for CTH. By testing for sensitivity to variations in strain- 
rate it will be possible to determine whether high strain rates change the final crater dimensions 
(ie., contribute to supralinearity as suggested but never proven by many investigators), or merely 
alter die dynamic shape and rate of development (as suspected by POD). 

The issue of high impact speed vapor momentum enhancement needs to be studied more 
accurately. POD's present analysis suggests only a relatively minor enhancement for a vaporizing 
projectile. Analysis by others (Lawrence, 1989) suggest a much larger effect at very high impact 
speeds where substantial portions of die target also vaporize. The issue is whether the larger 
momentum translates into noticeably larger craters (or a different shape) for a given impact speed. 

For additional work in the future, POD recommends that die analytic approach for transition 
between pure cratering (in semi-infinite targets) and marginal perforation be more fully 
developed. Additionally, an approach needs to be developed to determine the back surface 
spallation and die perforation hole sizes. Also, perforation studies need to be done for other 
materials. 

CTH calculations could be done to track out the data observed by HOrz of varying hole size 
versus target thickness. This would allow a firmer understanding of this behavior. HOrz's work 
is important because it indicates one of the very few techniques for deciphering perforation data 
for projectiles which are not much larger than the target thickness. 

Based on the ricochet and oblique impact laws reported here, studies should be done to extend 
die laws and correlate oblique impact crater asymmetries versus diameter and depth to allow 
direct interpretations of impact angles. 


64 



Studies should also be done for layered targets, such as die thermal blankets on LDEF, and th 
thermal paints on aluminum. 

Additionally, those individuals working on alternative scaling laws (e g , McDonnell, Cour-Palais 
should use die "pi" terms identified in this series of reports in order to provide better physics 
based fits to die data. 
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APPENDIX a - FIGURES 
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Figure 1-P:u Logic 1-D Verses Bernoulli Flow 
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Figure 2-P:u Logic Rigid Projectile Into Soft Target 
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Figure 3. Plot of the Petch Law’s Supralinearity 
Downgrading Factor vs. Crater Radius (for Aluminum) 
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Figure 5. Plot Showing How Watts’ Eq. 59 Varies 
With N (Al into Al Impact) 



Figure 6. Comparison of Watts’ Equations with 
Horz Data for Aluminum 6061 -T6 Impacts at 6 km/s 
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Figure 7. Plot Showing How Watts’ Eq. 59 Varie 
With N (Al into TFE Teflon Impacts) 
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Figure 8. Variations in Equations (82), (83) and (84} 
with Angle of Incidence - Theta 
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Figure 9. Variations in dc/dp 
with Angle of Incidence 
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Figure 10. Ricochet Angles for Al into Al 
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Figure 11. Plot of CTH Results Showing How P/dp Varies 
With Velocity and Projectile Density for Al into Aj 
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Figure 12. Plot of CTH Results Showing How dc/dp Varies 
With Velocity and Projectile Density for Al into Al 



Velocity (km/s) 


iiHNltlil IIII4I" iliilMWIII'l li l! |l I ffljif HMiittl ! ^b | l-«f lj i ilHl | -i!l> |- jfHi-HPH M PHi ; ^-H i ^ .HlM ! || i - iq H MM»j- silM ijblillWillM II II • IfiWl-Hil} |..| , 1. t fa* « .lip ll-il h(F ** M4I» » «• m-m + 


Figure 13. Plot of CTH Results Showing How P/dp Varies 
With Velocity and Target Yield Strength for Al into Al 
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Figure 14. Plot of CTH Results Showing How dc/dp Varies 
With Velocity and Target Yield Strength for Al into Al 
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Figure 15. Plot of CTH Results Showing How P/dp Varies 
With Projectile Density and Velocity for Al into Al 
10 


—Velocity 

-Velocity 

—Velocity 

—Velocity 


1 km/s 
5 km/s 
10 km/s 
15 km/s 


Density (g/cc) 


Figure 16. Plot of CTH Results Showing How dc/dp Varies 
With Projectile Density and Velocity for Al into Al 
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Figure 17. Plot of CTH Results Showing How P/dp Varies 
With Target Yield Strength and Projectile Velocity 
For Al into Al 
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Figure 18. Plot of CTH Results Showing How dc/dp Varies 
With Target Yield Strength and Projectile Velocity 
For Al into Al 
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Figure 22. 
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Figure 23. 

Al Projectile Into Aluminum ot 10 km/i 
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Figure 25 . 
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Figure 27. 
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Figure 30. 
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Figure 31. 
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Figure 33. CTH Data Showing How P/dp Varies With 
Projectile Density and Velocity for Al into TFE Teflon 



Figure 34. CTH Data Showing How dc/dp Varies With 
Projectile Density and Velocity for Al into TFE Teflon 
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Figure 35. CTH Data Showing How P/dp Varies Wit 
Target Yield Strength and Projectile Velocity 
For Al into TFE Teflon 
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Figure 36. CTH Data Showing How dc/dp Varies With 
Target Yield Strength and Projectile Velocity 
For Al into TFE Teflon 
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Figure 37. CTH Data Showing How P/dp Varies With 
Projectile Density and Velocity for Al into TFE Teflon 
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Figure 38. CTH Data Showing How dc/dp Varies With 
Projectile Density and Velocity for Al into TFE Teflon 
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Figure 39. Plot of CTH Results Showing How P/dp Vanes 
With Target Yield Strength and Projectile Velocity 
For Al into TFE Teflon 


-Velocjty 

Velocity 

Velocity 

Velocity 


Yield Strength (kbars) 


Figure 40. CTH Data Showing How dc/dp Varies With 
Target Yield Strength and Projectile Velocity 
For Al into TFE Teflon 
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Figure 41. CTH Data Showing How P/dc Varies With 
Projectile Density and Velocity for Al into TFE Teflon 
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Figure 42. CTH Data Showing How P/dc Varies With 
Target Yield Strength and Projectile Velocity 
For Al into TFF Teflnn 
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Figure 48. 
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Figure 49. 
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Figure 52. Comparison of Watts’ Equations 
with CTH Predictions for Aluminum on TFE Teflon 
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Figure 60. Comparison of CTH Results versus 
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Figure 92. Comparison of CTH Results versus 
Watts’ Eq. 59 for Al into TFE Teflon 
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Figure 124. 
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Figure 125 . 
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Figure 126 . 
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Figure 128. 
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Figure 130. Comparison of Watts’ Equations with 
CTH Predictions for Aluminum on Aluminum (Al 6061 -T6) 
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Figure 131. Comparison of Watts’ and 
Bruce’s Equations for Aluminum on Aluminum 
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Figure 135. Comparison of Watts’ and 
Dunn’s Equations for Aluminum on Aluminum 



0 1 2 3 4 5 6 7 


Velocity (km/s) 


A-71 





Figure 136. Comparison of Watts’ and 
Goodman-Liles’ Equations for Aluminum on Alurr^ium 
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Figure 137. Comparison of Watts’ and 
Herrmann-Jones’ Equations for Aluminum on Alumfium 
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Figure 138. Comparison of Watts’ and 
Sawle’s Equations for Aluminum on Aluminum 
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Figure 139. Comparison of Watts’ and 
Sedgwick’s Equations for Aluminum on Aluminum 
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Figure 140. Comparison of Watts* and 
Sorenson’s Equations for Aluminum on Aluminum 
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Figure 141. Comparison of Watts’ and 
Summers-Charters’ Equations for Aluminum on Aluminum 





Figure 142. Comparison of Watts and 
Bruce’s Equations for Aluminum on TFE Teflon 
^Assumes St Equals Yt f or TFE Teflon ) 
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Figure 143. Comparison of Watts’ and 
Cour-Palais’ Equations for Aluminum on TFE Teflon 



A-75 






Figure 144. Comparison of Watts' and 
Dunn’s Equations for Aluminum on TFE Teflog 
(Assumes Sigma-yt Equals Yt for TFE Teflon! 
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Figure 145. Comparison of Watts’ and 
Sedgwick’s Equations for Aluminum on TFE Teflin 
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Figure 146. Comparison of Watts’ and 
Sorenson’s Equations for Aluminum on TFE Teflon 
(Assumes St Equals Yt for TFE Teflon) 
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Figure 147. Comparison of Watts’ and McDoniell 
& Sullivan’s Equations for Al into Al (Al 6061 -T|) 
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Figure 148. Comparison of Watts’ Equation 59 versus; 
Other Investigators’ Equations for Al into Al 
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Figure 149. Comparison of Watts’ and McDonnell 
& Sullivan’s Equations for Al Into TFE Teflon 



Figure 150. Comparison of Watts’ Equation 59 versus 
Other Investigators’ Equations for Al into TFE Teflon 
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Figure 152. 


Soda Lime Glass->AIuminum 6061 -T6 
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APPENDIX 


Penetration Depth Equations 

Reference No. 27 : 

p/d = 2.28 (P D /P t ) 2/3 (V/C) 2/3 , 


V < 9 km /sec 


Reference No. 28: 


p/d = 1.96 (P p /P t ) 1/2 (V/C) 


2/3 


V < 6 km /sec 


Reference No. 29: 


p/d = 1.5 (P p /P t ) 1/3 (p p V 2 /2S,) 1/3 


V < 8 km /sec 


Reference No. 30 : 

p/d = 2.35 (p p /P t )°’ 7 ° (V/C) 2/3 


V < 9 km /sec 


Reference No. 31: 


9 1/3 

p/d = 0.63 (PpV^/ Oy^) i/ 


V < 7 km /sec 


Reference No. 32: 


.0.537 0.576 


p/d = 0.482 (p p /p t > (V/O 


(Y t /p t C 2 )' 0 ' 235 V < 21 km/sec 


Reference No. 33: 


p/d = 8.355 x 10 ’ p p t (V /H t ) 


~4 _ 2/3 _ ~l/3 
P 


V < 9.5 km/sec 


Reference No. 34: 


p/d = 2.0 (P p /P t ) 4 ' 52 (V/C) 


1.136 


V < 9 km /sec 
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(P-1) 


(P-2) 


(P-3) 


(P-4) 


(P-5) 


(P-6) 


(P-7) 


(P-8) 



Reference No . 35 : 


0.17 , it2 /c .0.285 


p/d = 0.311 (P p /P t ) * <P p V /s t> 


Reference No. 36: 


2/3 /. tt2/-q -v 1/3 


p/d = 0.36 <C p /P t > ' <P t V*/B,)*" . 


V < 6 ign/sec (F 


Reference No. 37 : 

p = 2.973 x IQ’ 7 d 1 * 1 H t '°‘ 25 p °* 5 P t '°‘ 167 V 4/3 , V < *3 Jcm/sec (F 


= 1.129 x 10 d 


u "t K P t 

-6 .1.056 u -0.25 0.5 -0.167 E -0.33 y 4/3 


Pp p t 


Crater Diameter Equations 


V < 'v 3 km /sec (F 


Reference No. 18: 


d 2 p/d 3 = 34 (p„/p t ) 3/2 (V/C) 2 , 


V < 4 kin /sec 


Reference No . 35 : 


.0.5 , XT 2/o v 0.845 


a d^p/d 3 = 0.120 Cp p /p t )“- : ' <Pp V ' S t > 


V < 7 -m/sec 


Reference No. 28: 


d h 2 p/d 3 = 30.25 <p /p,) 3/2 (V/C) 2 


V < 6 ;jn/sec' 


Reference No. 30: 


d 2 p/d 3 = 44.10 (P p /P t ) 2/3 (V/C) 2 . 


\ 7 < 9 *m /sec 


Reference No. 33: 


d 2 p/d 3 = 2.65 x 10‘ 9 p d 7/6 p,' 1/2 V 2 /H t . V < 9 5 km /sec 



Reference No. 36: 


a d h 2 p/d 3 = 0.16 <P p /e t )3/2 P p y2/B t ' V < 6 km/sec 

Notation 

d. ... crater surface diameter (cm) 
n 

d ... projectile diameter (cm) 

p ... crater depth (cm) 

2 

B t ... target material Brinell Hardness (dynes/cm ) 

C ... speed of sound in target material (cm/sec) 

E t ... target material elastic modulus (GPa) 

2 

H t ... target material Brinell Hardness Number (kg/mm ) 

2 

S ... target material static shear strength (dynes /cm ) 

2 

S t ... target material dynamic hardness (dynes/cm ) 

2 

... target material dynamic shear strength (dynes /cm ) 

V ... projectile impact velocity 

a ... crater shape factor 

a = 0.75 if p > d h /2 
a = 1.00 if p < d^/2 

3 

p ... projectile material mass density (gm/cm ) 

3 

p ... target material mass density (gm/cm ) 

' 2 
o ... target material dynamic yield strength (dynes/cm ) 

Material Properties 
10 2 

B t = 1.27 x 10 dynes /cm 

C = 5.10 x 10 5 cm/sec 
E t = 7.38 x 10 10 N/m 2 
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(C-6) 



S = 2.83 x 10 9 dynes /cm 2 
10 2 

S t = 6.37 x 10 dynes /cm 
9 2 

Y t = 2.78 x 10 dynes /cm 
3 

p =2.71 em/cm 
P 

3 

P t = 2.84 gm/cm 

= 1.85 x 10 10 dynes/cm 2 
H t = 130 kg/mm 2 
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HYPER VELOCITY IMPACTS ON SPACE DETECTORS: 
DECODING THE PROJECTILE PARAMETERS 
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Canterbury, Kent CT2 7NR 
United Kingdom 


ABSTRACT. Hypervelocity impacts in space have been used as a tool for the study of the paniculate environment 
throughout the space age. Detectors, both designed and 'incidental', have utilised cratering, penetration, 
momentum and (transient) plasma to detect such projectiles. Decoding the impacting projectile parameters from 
such target behaviour has often been limited by the calibration data available in the mass- velocity plane, and by the 
need for data to extend over a range of projectile dimensions from sub-micrometre to centimetre scales. LDEFs 
return to Earth with a wide variety of target materials and a very high definition of the particulate flux, along with 
its angular dependence, has provided the opportunity to accurately assess environmental impact data. Yet it forces 
the issue of which formulae are appropriate in this size range and which formulae can reliably extend over the wide 
range of velocities and particulate size regimes. The analysis of existing hypervelocity impact data for iron 
projectiles leads to ballistic limit penetration formulae which extend over orders of magnitude of target material 
densities and relative strengths as well as velocity. Although not explicitly calibrated for differing projectile 
densities, the form of the relationship (and established by other data) readily lends itself to the incorporation of the 
effects of projectile density. 


1- The Hypervelocity Impact Data 

Hypervelocity penetration studies on thin foils, of thicknesses between 0.8 pm and 4.8 um, have 
been carried out by McDonnell (1970, 1979) using a 2 MV Van der Graaff Accelerator. Since 
experiments were limited to using an iron dust source as projectiles (due to the nature of the 
electrostatic accelerator and dust availability) a range of impactor-target density ratios was 
investigated by varying the target material. The density range of materials explored ranges ftt)m 
mylar (a plastic polymer of density 1.395 gem* 3 ), to platinum foils of density 21.45 gem* 3 . The 
full data set, incorporating data by McDonnell (1970, 1979) investigates iron impacts onto various 
metallic and mylar foils between velocities of 1.0 kms* 1 to 16 kms* 1 . 

The experimental program is given by McDonnell (1970). For each impact event the particle 
diameter, d, the perpendicular impact velocity, V, the penetration crater hole diameter, Dh, and the 
foil thickness, f, were measured. Each of these impact events, on a (f/d)-V plot, is then associated 
with its penetration diameter, also normalised to the foil thickness, (Dn/f). For each projectile and 
wrget combination, the set of data obtained this way can then be plotted to define contours of 

(Dh/0- 

. P° r a given impact velocity and foil thickness, the crater diameter increases with the panicle 
size; As the particle size is. decreased, in the experiment the hole diameter too decreases, until a 
minimum detectable hole size is reached. This then represents the marginal penetration cut-off; 
the ballistic limit is defined by the contour approaching the asymptote of Dh/T » 0. These 
marginal penetration limits (very close to the ballistic limit) established for each projectile-target 
scenario are shown in Figure 1 as heavy lines. In some cases, the limiting threshold of the 
experimental technique was insufficient to define well this marginal cut-off, so that an upper limit 
of perforation could only be achieved. This is seen to occur where the smallest (Dn/f) ratio 
observed is not a heavy line in Figure 1. 

The equations of these penetration contours can be defined, and so, for each target and projectile 
scenario, a penetration limit established. Table 1 gives the results of this for each data set. The 
marginal perforation limits (f/d) are given as functions of velocity; the velocity regimes in which 
these are measured are also stated, along with their (Dh/ 0 limiting values. 




Table 1. The marginal penetration limits of (f/d) along with their (Dg/O limits ire define^or iron projectiles 
impacting the ten target materials between the velocities giveo. Other data used in this pSer for Aluminium 
targets extends to 16 km*-*. 1 


Projectile 

Material 

Iron 

Iron 

Iron 

Iron 

Iron 

Iron 

Iron 

Iron 

Iron 

Iron 


Target 
Material 

Aluminium 
Copper 
Stainless Steel 
Iron 
Silver 
Platinum 
Gold 
Titanium 
BerylliumCopper 
Mylar 


Marginal 
Perforation Limit 

f/d* 1.298V 0 - 586 

f/d=0.51 8V0-6W 

f/d*0.419V°- 753 

fld*0.517V°- 581 

I/d=0.599 VO- 708 

f/d»0.270V 0 -549 

£fd=0.542V°- 548 

f/d=0.58lV°- 737 

f/d«0.373V°- 8 9 4 

f/d=1.721 V 0 - 582 


ploctty Range 
1 (kms* 1 ) 
r i.6-4.4 

3 3.5-S.6 
- 2.S-3.7 
: 3.0-4.0 
! 4.2-5.3 
! 2.0-3.6 
11.2-3.5 
] 3.S-6.7 
15.1-6.5 
13.4-4.8 



velocity (knv%) 


Figure 1. Perforation limits for all target materials with their normalised Dfj/f values for irc~ projectiles 


2. Development of a Penetration Formula 

It is seen from the data in Figure 1 that the (Dh/ 0 contours are parallel ne^ to marginal 
perforation: noting the experimental measurements of aluminium (which extend injhis fashion to 
16 Jems* 1 ) and the similar behaviour of all these ductile targets over the limited ran A compared in 
Figure 1, it is likely that all these marginal perforation limits can be cxirapolat^ to very high 
velocities. A general equation describing all these proj'ec tile-target scenarios is, Jowever, first 
established in the measured range by incorporating certain properties of the projeefle and target 
materials. Looking at Table 1 and Figure 1 it is clearly seen that the (f/d) margipl penetration 
limit of the iron - * mylar combination is very different to the iron— iplatiman da ta^ This may be 
explained by the fact that platinum has a density some fifteen times higher than thayof mylar (see 
Table 2). This introduces the idea that a generalising equation of the complete dat Act must have 
a projectile-target density ratio scaling factor contained within it. Examining data ij|Figure 1 and 
the material properties given in Table 2, it is clear that there is a correlation between the (f/d) 
L - limiting contour and the target density, with the material density decreasing as (f/d) values 
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i - y-uncirv target material would reejuire a 
increase This is seen 10 be realistic, because a lowe ^ t0 a high density material, with 

" u,eri * 1 dtMiiy " ,d 10 “ inOT “ e or 

marginal penetration limit (f/d). f ctor a i onc would not completely satisfy the data. 

But it is also evident that a density scaling l0 ^ c i|- relative levels on Figure 1, and 

This is demonstrated by ordering lh * . ns ities Beginning with the greatest (f/d) values, the 

comparing this with *«. ° V ?^ (1 395 gcm 3 ). alumfnium (2.71 gP»\«“i“™£!*« 

materials are. with then (8.9Z8.2 gcm^). gold/iron/staii 



too ’high' relative to titanium. 

get materials. 

h Au 


Obviously, another material propert^neces^^^n Cy on material tensile strength, 

introduced into this ballistic htmt seen ^ matwial would imply a lower (f/d) line than 

On Figure 1. it can be envisaged that a high wength ^ strongcr the target is. the thinner 

that of a lower strength mat * n i° f * *^ atcr expansion ceases in proximity to the rear face, 
the target can be before the late sta g . end D f this marginal perforation process. If 

Spallation (involving renri/e available materials it is apparent that iron and 

one now considers the tensile *£*"** ., strengths along with titanium, whereas copper, 

stainless steel have rela^ely high ' “3£ff£*£^a»w*ito. Overall, if the targe, tens* 
beryllium copper, silver and gold have retaoveiyww^K^ fonnu j a d 0 n. then this appeus to explain 
strength is included into the general m Therefore, for a panicle of diameter, d 

the discrepancies that a simple density f p [’^-3], tensile strength, o T 

[cm], and density. pp[gcm*3], iropacn^ a target : of de i 5 g v [ kms -l], the marginal 

[MPa], and thickness, f [cm], at ^^^“^JS^SBneterl /l. B, C. and D are to 
perforation limit of the form given below “ prop ® £ „ that of aluminium, where 0*1 * 80 

fee determined. The target tensile strength is reterenceo 

MPa. 


l. A (— )*(— ) C v D 

d p T r 


(eqn 1) 


The form of the parameter D for the . cx P c P r !^?g ^^^Sc^ity^xponents are observed to 
Studying the marginal penetration c( l uatl °" s K^Tfmav be concluded that there appears to be no 
rangefrom between 0.S48 through to 0; 8 ^Vi^ s ^„a«d here. Furthermore. Figure 
correlation between these and their target ma P ' P*j Den etradon slopes (being the velocity 

1 demonstrates the^simiimtyteweenth It exp0 „ents have no significandy arg 



TOs U Sen die value. of * '”^„ti^s 7. dimugh ,0 

equation 1 iorri.i.and C » ^Ig^jSu !l“lvt Mth“ su'ge^xperi^mri bus due 
10 ), that this exponent and those in ^ at > . . „ ^^lin the measurement range, 

to the lack of a weak, but significant, dimensional seal ng wwn ^.^ation u uscd t0 obtain the 

Analysis using a three imensionalme a ^ n 1. This is calculated at a velocity o 

,st possible fit between the data > » Tibtel wM ch the measurements were taken. The 

fc rieS **»•* «- 


best po 
kms-i, : 
solution 
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'.0.635(^r‘(^) 0,: ‘v‘ 
d P T °T 


(cqn 2) 


, . , ■ ninned against the (f/d) date of Tab® 1 in Figure 2 at 

A plot of (f/d) as given by equation 2 » P ,k e B ecn eralised equation and Se data would sho\ 
velocity of 4 kins • A perfect m, 5* *“££« 2 stows lit « the genetaliseSquanoo compere 
all the data points on the straight line. Fig source data for the prc#crile-targei densn 

favourably with the data when one ■““*£ <£ a SmStSe; soSrcSata for the densir 
ratios and tensile sjengths ^ cove ^ strength, compared to theleferenced value c 

ratio ranges from 0.37 »o 5.64. wnus the fi| was calculated at th«nean experiments 

aluminium, ranges fromO.l 3to2^u. s describe the actual data.fus equation migh 

^om^cSroos, be used .. high* vcicc.de, j 


A 5 4 VA (eqnj 




SSrS^ #Ti 

F^CU^Ag 


jd (data] 



„ o, . « * « 2 3 3j 4 r 3 

Figuic2. Cofltparimo between Ike Cat* tmtequeoeo 2. | 

The foil IhieknKses used 

i nsw! ^ssSfSssss s«cr*:5Ki= & 

equation 2, desenbes the scaling factor, validating #s equation up *i 

extended to incorporate a P«jc Unifying the microscopic and microscopic margina 

millimetre sized projectiles and th^fore wjuyi^g h clocitj r impact dsf at both these size 
penetration regimes. This is achie a y P . winarion 2 anolics to Ditnecules with a mcai 

ranges The McDonnell dau pra semed tn ^ potion depih (P) fo 

pla^tMoTe £ojSd£ d2u2w. for projectiles of 5 mm in di.nte.er, yifmg: 

P .„/P r„n.5mm, teqo 3J 


0 0.667 y 0.667 

£- M5(£E) <i> 

d Pr 


If these two equudons ■"Tmw^tM^e^rfl^nd^^imi^twe 

leads to f/d - 2.648 for d - 2 1® ^. .. referring to opticSy detf table perforations 
C, is taken as 5.105 kms*. The ^ -knee" in the hole growth curvejs ^c parocle sm 

corresponds for aluminium targ . of thc f 0 y thickness, whereas thjstnct balbsoc to 

increases, namely a value ^ may thereby result, but thejcomparison cleari 

l&^d^SSfSS — .h. .« * * r* r °* over fc * 
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orders of magnitude in panicle size. From this comparison a scaling factor index is calculated by 
interpolating linearly on the logarithmic scale: 

index = [ log (3.894 / 2.648) ] / [log (5000 / 2) ] = 0.049 (eqn 4) 

This value would represent the exponent of the particle dimension, taking on the form: 


f/d«d 


0.049 


(eqn 5) 


(all other parameters constant) 

This is seen to compare very favourably with the dimensional scaling factor given by Fish and 
Summers (equation 8). and used by Naumann (equation 10). and Cour-Palais (equation 1 1), 
namely 0.056. The Pailer and Grfln equation (9) by contrast gives a value of 0.2 (see next secnon). 
The simple dimensional scaling we calculate shows such close agreement with a weUesu Wished 
and widely used -value, that we must adhere to this accepted value of 0.056. The resultant 
divergence between indices of 0.049 and 0.056 over 3 magnitudes in panicle diameter is less than 
5%. Applying this factor to the penetration formula of equation 2 expressing d in cm and 
normalising for this term, gives a relationship of the form shown in equation 6. where the panicle 
diameter is scaled to 2 pm: 


1.023 d 
d 


.0.056 


o 0.476 a 0.134 „ _ 

Pp \ f °A1 \ , .0.664 


( 2 ) ) 


Pt 


(eqn 6) 


In a further, but very important, consideration we must note that although equation 6 would 
represent and scale the microparticle data acquired from pm to cm dimension data, it is never the 
less acquired from a finite range of panicle diameters. The electrostatic accelerator projectile 
"mass spectrum" which is used for this data set, for example, shows dependency on velocity, 
namely m « 1 / V 4 and hence because d •« m 1/3 we find that the panicle diameter •* 1 / V 4 «. 
What is important, however, is the change in panicle sire along the observed Cfd contour. The f/d 
values at low velocities refer to measurements from larger particles than at higher velocities. This 
change of scale for a set of measurements on one foil is, therefore, dictated by the measured 
contour itself (namely by the particles selected) rather than by the spectrum of available particles. 
This effect of dimensional scaling within the measurement range thus leads to a bias to any of the 
observed contours of marginal or supra-marginal perforation, and leads to an apparent reduction in 
the velocity exponent The apparent dimensionally independent velocity exponent (i.e. of the f/d 
locus for constant particle sire d) can be shown to be 1/1.056 of its 'true' value. The 'true' value 
average value for this set of 10 target materials surveyed is thus 0.664 x 1.056 « 0.701 in a 
dimensionally scaled formulation. We must re-norroalise (6) with this exponent change and finally 
we have: 


-f« 0.970 d 0 056 
0 



(eqn 7) 


This equation (7) is termed McDonnell 1992A and is the scale corrected version of equation 
6 (Sullivan 1992); Vis in km s* 1 and d in cm. The results arc compared to other penetration 
relationships in Figure 3; the equation yields very favourable agreement with experimentally 
determined ballistic limits for particles of varying dimensions, and velocity;!! also encompasses a 
wide range of target densities and tensile strengths and in decoding likely projectile parameters 
from meteoroids and space debris is not likely to be in error by more than some 10% for arbitrarily 
chosen parameters. 

For specific projectile target combinations where calibration data is available at appropriate 
velocity, the source data must be used. For example calibration data specific to iron projectiles 
impacting aluminium microfoils for Dh “fat velocities up to 16 km s* yields f/d ** 0.79 V • ® 
(McDonnell 1969). We can now update this with better insight into the dimensional scaling 
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bias within this data which will permit applications to dimensions outsype the range of 
measurement. This leads to dimensionally scaled formulation normalise^ to the Itikms* 1 * 
measurement which tracks better the velocity dependence for this particular targjj configuration: 


f/d * 1.272 V 0 - 806 d 0 - 056 I (eqn 8) 

This equation is termed McDonnell 1992B and agrees well with thj experimentallj 
determined ballistic limits of f/d ■ 6.55 at 16 km s’ 1 for Dh /f *1 and d=0.24pm ^ _ 

If we compare the results of equation (7) for iron panicles impacting aluminium to the 
specific iron onto aluminium formula of equation (8), we find; 

f/d = 0.970 d 0 056 (7.8/2.7) 0 - 476 V 0 - 701 

« 1.607 d 0 05 ^ V 0 - 701 « 1 1.22 d 0 056 at V= 16km s* 1 i (eqn 9) 

c.f. from equation (8)- f/d - 1.272 d 0056 V 0 - 806 m 1 1.88 d 0 - 056 at V= I5kms | 

We note a modest divergence, from the application of the generalised eq^tion; for other 
conditions this divergence could be greater, and recommend, where available, fjrmulae derived 
from data closest to the panicular impact conditions. In further developmtyj of penetration 
ralationships, and taking note of the strength and density functional rcladonshipfifound in the 10 
target survey, we can extend the iron onto aluminium data, which now incorporjjes dimensional 
scaling, to cope with different target strengths and densities and (but with lessficcuracy in the 
lattter ) for different projectile densities. This yields : • 





(eqn 10) 


3. Previous Penetration Formulae; Comparison. 

Other widely used equations are maybe usefully compared, having been convened to the same 
units as used to define equation 1 and with d in cm. The foil ductility, e, is dimt|sionless, and a 
is the impact angle between the target normal and projectile trajectory. We compaf^ for example 

1 - 0.79 V 0,763 McDonnell (*79) (eqn 1 1) 

G - 


^ ■ 0. 57 d 0-054 e°' 03< ( — )°* 3 

d Pt 


i-0.772d 0 V'%JV' ! 


v a«75 
(V cosa) 0,M 


Fish & Summers (*65) (eqn 12) 


Pailer & GrOn :*80) (eqn 13) 


f 

d 


.0.056 

a 


0.32 v 0.«75 

Pp v 


Naumann (*66) (eqn 14) 


» 0.635 d°‘° 56 pp V Cdur-Palais (if 79) (eqn 15) 

The McDonnell (1979) equation (as described) was used, because of its validaqjn at velocities 
higher than light gas gun data, to extrapolate to 68 kms* 1 for application to the Gjjtto probe with 
its comet P/Halley encounter in 1986. Because this equation explicitly described Jon projectiles 
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»• _ r, r tnr was used; but this has 

^■-rr^.-sw 

H^v CoSsons can now be made between the^ pwa^n d ^ cnsiona i bias agree 
SSnentoMhe equations developed when now cotrected Thc paniclc and target 

4 to* . soled ta iron projecules JSton Eon (SMM) Ooonto *nd 

and Griin equation, used to interpret data from tn nerforarion equations. It can also be seen 

Brownlee! 1986), does not generally agree with Opi^ ^ ^dimensionally scaled equation 
that overall equations by Fish and Summers, Cour • ^ ^ velocity and projectile size 

- - - ih * d “ 1 

scaling now introduced. 


Scaled for Fe > Al at 4km/s for 2 micron panicle 



N (eqn 14) 

| P AG (eqn 13) 

. - - C-P (eqn 15) 

„/ McD A S (eqn 6) 

' McD 1992A (eqn 7) 

• s / McD 1992B (eqn 8) 
'McD l992C(eqn 10) 

FAS (eqn 12) 


0.0001 


0.001 0.01 

Projectile diameter (cm) The eauation developed here 

Pi,„,. Marginal perforation equation, as a function of pcoj^nlc diameter, m*. 

is shown as the dashed line. 


P~-SS- 3 S« 

sample of thc McDonnell Uvin fo ^" e ^°° c ^ io n Equation depends on the closeness of an 

SSSS^Si SSfSafiUS. d„. evaiUble. Table 3 shows values calcula.ed for 
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different formulae to two situations where impact calibration is available at# and 16 kms' 1 . N 
of guidance for the choice of appropriate equations from the new formulations are also given. 


ues tor he onto Al 




X= 1.023 d 0 056 (jk)°' <76 (SilJ 0134 v°-< 


X = 0.970 d 0056 f-^l ° 476 


yO.701 


1=0.79 V 0 - 763 
d 


X = 1.272 d 0 056 ^* 06 
d 


corn’d.. (5Al)° 134 V 0 M 6 




Table 3. Comparison of penetration formulae developed in this work. Values used for jpnputaion area:p yy * 
2.71 (gfcm 3 ). pFc ■ 7.87 (g/cm 3 ); d On cm) yields; (2pm)0°S6 ■ 0.6207; (0J4jim) 0 ^ « 0 .5512; 
(PF^PAJ) 0 - 47 * -1-661 | 

McDonnell 1979 feqn ll). Iron onto Aluminium; no dimensional scaling iM velocity exponen 
biased by dimensional range in data; calibrated at 16 kms* 1 for dp =0.24 pm. | 

McDonnell 1992B feqn 8). Extends McDonnell 1979 to remove scaling |>ias in data and ti 
dimensionally scale. Applies therefore to a wide range of dimensions (e.g. Imogens to centimetres 
and velocities from 4 to 16kms‘* for iron projectiles onto aluminium. 

McDonnell & Sullivan 1992 feon 6) Covers wide range of metallic target strgtgths and densite; 
(including mylar) but has residual errors for some materials calibrated; velodcj range is ^kmr 1 
note that the velocity exponent is an average for the 10 targets and is biast| by the dimension 
range within the data. 

McDonnell 1992 A feqn 7). As McDonnell & Sullivan 1992, but removes dimensional bias 
within the measurement range. Applies to a wide range of target strengths Mjl densities, though 
the velocity exponent is the 10 target average and also small residual errors for jome materials will 
exist; therefore applies to wide range of velocities (4 kms* 1 to 16 kms- 1 at least) and dimensions 
(microns to centimetres). 


targets but includes removal 


Derived from McDonnell 1979, for iron prqjcctujs onto aluminium 
of dimensional bias, the inclusion of dimensional scaling and the 



















functional fonn for v f ation ofpan.clc 

i?o^c o7po™™y ‘“scale to arbitrary projecrile-targe. configurations which may 
"^a^'a.^ 

and quannfied and yet is unc aP la ' nc d. , {,iclt depends upon strain rate. Our functional 

accounting £3*2 ™X£££) stre^h separately could therefore 

KdlSS farmwid a si depended, strength reUdonship. 
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